ADA034845 


028096 


NRL  Memorandim  Report  3408 


Vectorized  Tridiagonal  Solvers 


J.  P.  Boris 


Plusmti  Dynamics  Branch 
Plasma  Physics  Division 


This  work  was  sponsored  by  the  Office  of  Naval  Research 
under  Contracts  RR- 1 4-03-02  and  RR-1 1 -09-4 1. 


REPRODUCED  BY 

NATIONAL  TECHNICAL 

information  service 

U.S.  DEPARTMENT  OF  COMMERCE 
SPRINGFIELD,  VA.  22161 


NAVAL  RESEARCH  LABORATORY 
Wadiiaftom  D.C. 


Approved  for  public  release:  distribution  unlimited. 


security  classification  of  This  pace  Dmtt  Enfrmd)  


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I REPORT  number  i2.  aOVT  accession  no.]  3-  RECIPIENT'S  CATALOG  NUMBER 


REPORT  DOCUMENTATION  PAGE 


NRL  Memorandum  Report  3408 


4.  title  (and  Sublllld) 


VECTORIZED  TRIDIAGONAL  SOLVERS 


S.  TYPE  OF  REPORT  » PERIOD  COVERED 

Interim  report  on  a continuing 
NRL  problem. 


6.  PERFORMINO  ORG.  REPORT  HUMBER 


*.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


Naval  Research  Laboratory 
Washington,  D.C.  20375 


S.  contract  OR  GRANT  NUMBeRfF; 


10.  PROGRAM  element.  PROJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 

NRL  Problems  H02-39  & H02-51 
RR140302  and  RR110941 


II  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 
Arlington,  Virginia  22217 


IJ.  REPORT  DATE 

November  1976 


13.  NUMBER  OF  PAGES 


*■  MOMITORINC  AOeNCY  NAM6  ft  AOORESS<<f  aittfnt  from  ControUing  Ofticm)  SECURITY  CLASS.  r«popi; 

UNCLASSIHED 

ISi.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  fa/ Ifila  RFporo 

Approved  for  public  release;  distribution  unlimited. 


IT.  DISTRIBUTION  STATEMENT  (ol  aSalracI  mimnd  In  Block  20,  II  dllloconl  tram  Rmport) 


If.  SUPPLEMENTARY  NOTES 

This  report  was  sponsored  by  ONR,  Under  Project  RR140302,  title  Numerical  Hydrodynamics, 
and  Project  RR110941,  title  Computational  Fluid  Dynamics. 


19-  KEY  WORDS  fConNnti*  on  1/  Identify  by  ttoek  number) 

Tridiagonal  systems 
Parallel  Computation 
Vectorization 


20.  ABSTRACT  (Contlnu0  on  rovoroo  o/rfo  ff  noeoooorr  «nd  td0nt/fy  by  block  number} 

This  report  briefly  documents  four  related  tridiagonal  solvers  optimized  for  parallel  computa- 
tion on  the  ASC  computer  at  NRL.  Both  single  and  double  precision  versions  are  available.  The 
report  also  contains  listings  of  the  routines  and  a test  program  with  results  for  a standard  test 
problem. 


EDITION  OF  I NOV  El  II  OBIOLETE 
S/N  0102-0I4-E60I 


security  CLAIIIFICATION  of  This  page  f-Bhon  Dotm  Bnttrod) 


CONTENTS 

I.  INTRODUCTION  1 


CONTENTS 

I.  INTRODUCTION  1 

II.  SOLUTION  METHOD 2 

III.  USAGE 4 

IV.  TIMING  6 

REFERENCES  11 

APPENDIX  A — Four  Double  Precision  Tridiagonal 

Solvers  for  the  ASC  12 


APPENDIX  B — Test  Program  for  Double  Precision  Solvers 

Including  Both  Single  and  Double  Precision  Results  ...  24 


Preceding  page  blank 

iii 


“A-1  * »A  * =A+i  - '•i 

for  i = 1,  2,  3}  Here  is  the  number  of  columns  (and 

rows)  in  the  corresponding  matrix  representation.  These  appli- 
cations include  implicit  diffusion  equations,  splines,  Poisson 

Equation  solvers,  and  many  others.  References  include  Roache^ 

2 

auid  Richtmyer  and  Morton  . 

BoundaJT  conditions  for  the  system  depend  on  the  coefficients 

a^  and  c^^^  . Vlhen  a^  and  c^  are  zero,  the  first  and  last  equation 
c c 

couple  only  two  of  the  unknown  values  {X^}  and  the  usual  double 

sweep  algorithms  work.  The  periodic  system  in  which  h Xjj  and 

Xjj  = X^  also  occurs  often  in  real  applications.  Both  periodic 

c 

and  aperiodic  boundary  conditions  are  handled  as  if  periodic  in 
the  routines  described  below.  Distinction  between  these  two 
conditions  occurs  only  through  the  values  assigned  to  the  input 
coefficients  a^  anc  c^j  . 

Eight  different  user-called  routines  are  provided  and  docu- 
mented in  the  SPL  library,  four  for  single  precision  and  four 
for  double  precision.  Four  extra  routines  are  required  for 
internal  use  by  the  vectorized  tridiagonal  solvers;  these  four 
are  not  to  be  referenced  by  the  user.  Two  of  these  auxiliary 


Note:  Manuscript  submitted  October  29,  1976. 
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routines  are  for  the  double  precision  solvers,  and  the  other  two 
are  their  single  precision  equivalents.  Appendix  A lists  the 
four  double  precision  tri diagonal  s'^lvers  and  the  t\ro  auxiliary 
coefficient- folding  routines  used  by  the  two  vectorized  solvers 
TRIDDV  and  TRIDDM.  The  single  precision  equivalent  routines  are 
not  listed,  being  derived  from  those  in  Appendix  A simply  by 

1.  Chajiging  the  D to  an  S in  the  fifth  letter  of 
each  routine  name,  and 

2.  Removing  the  EXPLICIT  REAL*8  statement  from  each 
double  precision  routine. 

Appendix  B contains  the  SPL  double  precision  test  program  and 
its  output  for  verifying  the  results.  The  single  precision 
output  is  also  included  but  not  the  corresponding  single  precision 
test  program. 

All  of  the  tridiagonal  solvers  in  this  series  are  named 
according  to  the  following  conventions: 

1.  The  first  four  letters  are  TRIP. 

2.  The  fifth  letter  is  S or  D for  single  or  double 
precision. 

3.  The  sixth  letter  denotes  the  purpose  and  usage  of  the 
routine  as  detailed  in  Table  1. 


II . Solution  Method  | 

All  discussion  of  the  algorithms  used  will  be  based  on  the  1 

single  precision  solvers  although  it  holds  equally  well  for  the  3 
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double  precision  solvers  with  the  fifth  character  of  the  six- 
character  name  changed  from  S to  D.  The  basic  scalar  solver 
was  originally  written  by  C.  'I'Tagner  under  the  name  TRIDAG  and 
employs  a modified  double  sweep  technique  to  allow  the  non-zero 
corner  elements  of  a periodic  boundary  condition.  The  current 
routine  TRIDSS  (scalar)  is  based  on  TRIDAG. 

The  vectorized  algorithm  is  based  on  NRL  Memorandum  Report 
3l4b  in  detail  although  the  basic  ideas  have  been  around  for 
years.  ’ The  number  of  equations  is  reduced  roughly  by  a factor 
of  two  by  substituting  successive  even-numbered  equations  into  all 
of  the  odd-numbered  equations.  This  folding  is  repeated,  in  TRIDSV 
(vector),  over  and  over  a^ain  until  the  residual  tridiagonal  system 
is  so  short  that  the  TRIDSS  scalar  algorithm  is  the  most  efficient 
solution  method.  The  vector  operations  which  result  over  the 
entire  system  of  equations  more  than  make  up  for  the  extra  overhead 
of  the  folding  and  unfolding  operations. 

TRIDSV  also  manages  memory  so  that  the  operands  are  contiguous 
in  memory,  thus  majcimizing  memory  access  speed.  The  price  of  this 
is  the  necessity  for  every  vector  entered  to  TRIDSV  (or  TRIDSM)  to 
be  at  least  twice  as  large  as  needed  to  hold  the  data.  Thus  a 
vector  of  length  2N^  would  have  to  be  passed  to  the  vector  routines 
but  only  the  first  locations  are  filled  with  data.  The  upper 
half  of  every  array  is  used  as  scratch  and  the  original  contents 
are  destroyed. 

TRIDSR  (repeated)  is  a simple  rewrite  of  TRIDSS  in  which 
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seversLl  different  tridiagonal  systems  are  solved  in  parallel  by 

bringing  the  appropriate  DO  loops  inside  TRIDSS.  The  operations 

are  then  all  vectorized  and  contiguous  across  the  N tridiagonal 

s 

systems  although  the  scalar,  recursion  loops  of  TRIDSS  are  still 
in  evidence  as  outer  loops. 

TRIDSM  (multiple)  is  a vectorized  rewrite  of  TRIDSV  corres- 
ponding to  the  vectorized  TRIDSR  rewrite  of  TRIDSS.  In  TRIDSM, 
however,  both  inner  and  outer  loops  vectorize  so  it  is  the  most 
efficient  of  all  the  routines.  The  next  section,  which  describes 
the  actual  usage,  should  clarify  this  apparent  jungle  of  routines. 

In  all  of  the  routines  the  usual  divide  by  the  b^  coefficients 
is  retained  as  minimizing  the  number  of  operations  even  though  one 
of  the  operations  is  an  expensive  divide.  Therefore  the  diagonal 
terms  cannot  be  zero  (generally).  The  simple  (but  maybe  not  alwavs 
necessary  or  sufficient)  requirement  of  diagonal  dominance  should 
prevent  an  undesired  zero  divide.  This  choice  was  dictated  by 
simplicity  and  flexibility. 

III.  Usa^e 

Let  A,  B,  C,  D be  the  Fortran  arrays  in  which  the  coefficients 

stored  by  the  user,  and  X the  Fortran 

array  where  the  user  '/reints  his  results  {x^}  stored.  Let  SCA  and 

SCB  be  scratch  arrays  of  the  same  dimensionality  as  X.  Let  N be 

c 

the  ntmber  of  equations.  The  user  should  ensure  that  A,  B,  C,  D,  X, 
SCA,  and  SCB  all  have  dimensions  of  at  least  2*N  . The  vector 
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routines  use  (and  despoil)  the  upper  locations  in  each  array. 

A single  tridiagonaJ.  system  can  be  solved  by  the  statement 
CALL  TRIDSV(NC,A,B,C,D,X,SCA,SCB) 
or,  for  double  precision  arguments  throughout, 

CALL  TRIDDV( WC , A , B , C , D , X , SC A, SCB ) . 

TSIDSV  automatically  calls  TRIDDS  for  greater  efficiency  if  NC  is 
small  auid  hence  vectorization  is  inefficient.  Thus  the  user  gene- 
ra!Lly  need  not  touch  TRIDSS,  TRIDSF,  TRIDDS,  TRIDDF  at  all.  In  the 
event  that  the  user  cannot  afford  the  scratch  area  in  the  upper  half 
of  each  airray,  he  can  call  instead  TRIDSS(lTC,A,B,C,D,X,SCA,SCB) 
where  now  the  arrays  need  only  be  dimensioned  N^.  Figures  la  and  lb 
show  what  speed  the  user  loses  in  going  to  the  scalar  altorithm  to 
save  space. 

In  multidimensional  cases  let  A,  3,  C,  D,  X,  SCA,  and  SCB  be 
dimensioned  A(ND,  2*NC)  where  the  first  subscript  loops  over  the 
different  tridiagonal  systems,  the  lower  half  of  the  second  sub- 
script contains  the  data,  and  the  upper  half  of  the  second  index 
is  left  free  as  scratch  space.  This  system  of  multiple  tridiagonal 
systems  can  be  solved  by  the  statement 

CALL  TRIDSM(ITC,A,B,C,D,X,SCA,SCB,NS,ND) 
where  NS  is  the  number  of  systems  (out  of  ND)  which  are  actually 
solved.  For  double  precision  real  arguments  throughout  the 
statement  is 

CALL  TRIDDM(NC,A,B,C,D,X,SCA,SCB,NS,ND) 

TRIDSM  automatically  calls  TRIDSR  x^hen  the  residual  systems  become 


sufficiently  short  (NC  ~ 12).  Since  vectors  in  TRIDSM  are  two- 
dimensional,  however,  the  hreakeven  point  for  efficiency  "between 
TRIDSM  and  TRIDSR  is  at  lower  NC  than  it  is  for  TRIDSV  and  TRIDSS. 
"■/hen  NS,  the  mamber  of  systems,  is  very  large  (>  200),  the  fact 
that  the  scalcir  algorithm  in  TRIDSR  has  20^  fewer  operations  than 
the  vector  algorithm  will  begin  to  show  as  better  performance  for 
TRIDSR.  TRIDSR  should  also  be  used  when  the  factor  of  two  scratch 
space  required  by  TRIDSM  is  prohibitive. 

In  the  great  majority  of  cases,  however,  the  most  efficient 
method  is  to  use  TRIDSM.  In  general  the  user  need  only  remember 
TRIDSV  and  TRIDSM  and  their  nearly  identical  argument  lists  which 
differ  only  in  the  obvious  necessity  to  specify  the  second  index 
dimension  and  span  in  the  case  of  multiple  tridiagonal  systems. 

rv.  Timing 

Table  2 contains  approximate  formulae  for  the  execution  time  of 
the  six  most-likely- to-be-used  routines.  Figures  la  and  lb 
graph  the  same  data  for  the  single  and  the  double  precision  routines 
separately.  The  fastest  solvers  are  the  vector  multiple  tridiagonal 
solvers  TRIDSM  and  TRIDDM.  The  single  precision  version  is  the 
fastest  for  obvious  reasons.  There  ane  23  operations  needing  about 
30  "clocks"  of  computation  per  point.  This  is  about  2.4y;3ec  per 
point.  Since  there  are  two  arithmetic  pipelines,  however,  and 
many  of  the  vector  operations  can  be  executed  back  to  back  with 
others,  the  speed  is  substantially  faster  than  could  be  obtained 
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■td-th  one  pipe. 


Table  2,  in  addition  to  presenting  formulae  which  quantify 
Fig.  1,  contains  the  size  of  the  routines  in  words  (decimal). 
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Fig.  1 — Timing  data  for  optimized  tridiagonal  solvers,  (a)  the  three  principal 
single  precision  solvers  (TRIDSS,  TRIDSV,  TRIDSM).  (b)  the  three  principal 
double  precision  solvers  (TRIDDS,  TRIDDV,  TRIDDM). 
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Appendix  A 

Four  Double  Precision  Tridiagonal  Solvers  for  the  ASC 
The  routines  listed  directly  in  this  appendix  have  been  optimized 
for  and  run  on  the  ASC  system  at  URL.  Self-contained  documentation  appears 
in  comment  statements  to  these  routines.  To  change  to  single  precision, 
alter  the  subroutine  names, as  described  in  the  text  of  this  report  and 
remove  the  "IMPLICIT  3EAL*8"  card  found  at  the  beginning  of  each  routine. 


Note  that  all  intrinsic  double  precision  function  names  must  be  changed  as 
well. 
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SUBROUTINE  TRIOOS  (NC>  B>  C>  0,  X,  Rb,  AINVQ) 


TRIDOS  (NC,  a,  Bf  C,  0,  X,  RB,  AINVQ)  CLASS. 

ORIGINATOR  - C.  WAGNER  CODE  7750,  NRL  197« 

LATEST  MOOS  - J.  BORIS  CODE  7706,  NRL  NOV.  1976 

DESCRIPTION:  THIS  ROUTINE  SOLVES  IN  DOUBLE  PRECISION  THE  TRI- 
OIAGONAL  matrix  EQUATION  M*X  = 0 WHERE  A(I),  B(I),  AND  C(I)  ARE 
the  only  three  non  zero  ELEMENTS  OF  ROW  I.  SPECIFICALLY  WE  SEEK 
X(I)  IN  AtI)*XlI-l)  ♦ B(i;*X(l)  + CCI)*X(I+1)  = 0(1).  FOR  PERIODIC 
BOUNDARY  CONDITIONS  X(l)  i XCNCtl)  AND  X(0)  = X(NC).  THE  APERIODIC 
CASE  IS  HANDLED  BY  SETTING  ACl)  AND  C(NC)  TO  ZERO,  THE  SCALAR  AL* 
GORITHM  EMPLOYED  HERE  FOLLOWS  FROM  SEPARATING  M INTO  A PURELY 

trioiagonal  Part  ano  a vector  in  column  nc  and  one  at  row  nc. 


ARGUMENTS:  (n  ARR(... 

NC  integer 

A R ARR(NC) 

B R ARR(NC) 

C R ARR(NC) 

0 R ARR(NC) 

X R ARR(NC) 

RB  R ARRINC) 

AINVQ  R ARRtNC) 


) IS  USED  TO  DENOTE  REAL  ARRAY  I DIMENSIONS) 


0 EONS  IN  THE  TRIOIAGONAL  SYSTEM  I 
COEFFS  OF  X(I-l)  IN  I-TH  EQN.  I 
COEFFS  OF  X(I)  IN  I-TH  EQN.  I 
COEFFS  OF  X(I+1)  IN  I-TH  EQN,  I 
0(1)  - INHOMOGENEOUS  PART  OF  I-TH  EQN  I 
SOLUTION  VECTOR  SOUGHT  0 
SCRATCH  ARRAY  FROM  CALLING  PROGRAM  $ 
SCRATCH  ARRAY  FROM  CALLING  PROGRAM  S 


LANGUAGE  ANO  LIMITATIONS:  THIS  ROUTINE  IS  WRITTEN  IN  OOuBLE  PRE- 
CISION FORTRAN  FOR  THE  ASC  (69  SITS).  DIAGONAL  DOMINANCE  IS 
STRONGLY  RECOMMENDED  FOR  GOOD  ERROR  CONTROL.  THIS  ROUTINE  IS  REC> 
ommenoeo  when  short  (NC  < 15)  trioiagonal  systems  are  BEING'SOL- 
VED,  THE  VECTOR  TRIOIAGONAL  SOLVER  TRIDOV  CALLS  TRIOOS  TO  SOLVE 
THE  SHORT  RESIDUAL  SYSTEMS  GENERATED  BY  THE  VECTOR  FOLDING  OPER- 
ATIONS. THE  REPEATED  TRIOIAGONAL  SOLVER  TRIOOR  IS  BASED  ON  AN  EX- 
PANSION OP  TRIOOS, 


TRIOOS  has  no  additional  user  entries  and  calls  no  auxiliary  sub- 
routines. 


c 


IMPLICi;  REAL*8  (A-M,0-Z) 

OIMENSZON  A(NC},  B(NC),  C(NC)r  0(NC)/  X(NC)«  R8(NC}f  AINVQCNC) 
NM  s SC  - 1 

00  10  I 3 1.  NC 

IF  (BCD  .EQ.  0.0>0)  B(l)  s 1, 00*70 

continue 

RBU)  X 1.0^0/Bm 
00  15  I X 2,  NH 

RB(l)xl.0>0/  (B(I)  - A(I)  * C(I-l)  • RB(I-D) 

INVERT  TRIOIAGONAL  PART  ON  Q 
AINVOd)  3 A(l)  • RB(1) 

00  20  I X 2,  NM 

AINVQ(I)  X - A(I)  • AINVQCI*!)  • RB(I) 

AINVQ(NH)  X AINVQ(NH)  t CtNH)  * RS(RH1 

00  30  J X 2,  nM 

1 X NC  - J 

AINVOd}  3 AINVOCI}  - cm  * RB(I)  • AINVO(Itl) 

RCAQxl.O^O/  (B(NC)  • C(NC)  • AINVQ(l)  - A(NC)  * AINVO(NM)) 

ENO  PRECOMPUTATION.  THE  FOLLOWING  ARE  ALL  THE  CALCS  INVOLVING  0 
Xd)  X Od)  * RBd)  ' 

00  40  I 3 2,  NH 

X(I)  X (0(1)  - A(I)  * Xd-m  * RBd) 

00  50  J X 2,  NH 

1 X NC  • J 

Xd)  X X(I)  - C(I)  • RBd)  * Xd*l) 

X(NC)  X (0(NC)  - C(NC)  * Xd)  - A(NC)  * X(NH))  • RCAQ 

OB  60  I 3 1,  NH 

X(I)  X Xd)  - X(NC)  * AINVQ(I) 

RETURN 


entry  BOtCHO  (N,  M,  It,  Jt,  A,  B>  C>  0.  X) 
HANOLES  compiler  bug  in  TRIOOF,  TRIOOG. 
RETURN 
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CO 

CO 

CO 

CO 

CO 

CO 

CO 

c 


SUBROUTINE  TRIOOV  (NCr  B,  C,  D,  X,  SCA,  SCB) 


TRIOOV  (NC»  A,  B«  C«  Of  Xt  SCA,  SCB)  CLASS. 

ORICINATOR  • J.  BORIS  COOE  7706,  NRL  JULY  1976 

OESCRIPTIONi  this  routine  solves  in  oouble  precision  the  TRI- 
OIAGONAL  MATRIX  EQUATION  H*X  s 0 WHERE  A(I),  B(I),  ANO  C(I)  ARE 
THE  ONLY  THREE  NON  ZERO  ELEMENTS  OF  ROW  I.  SPECIFICALLY  WE  SEEK 
X(I)  IN  A(I)*X(I-1)  ♦ B(I)*X(I)  ♦ C(I)*X(I+n  X 0(1).  FOR  PERIOOIC 
BOUNDARY  conditions  X(1)  x X(NC«1)  ANO  X(0)  x X(NC).  THE  APERIODIC 

CASE  IS  handled  by  setting  A(1)  and  C(NC)  to  zero,  the  algorithm 
used  here  follows  from  repeated  cyclic  reduction  (NRL  MEMO* 

RANOUM  RErdWT  3104,  OCT.  76).  THE  AUXILIARY  ROUTINE  TRIOOF  IS  USED 
TO  AID  WITH  STORAGE  FOR  THE  NX  COMPLIER.  ALL  ARRAYS  MUST  HAVE  AT 
LEAST  2*NC  LOCATIONS  TO  ALLOW  FOR  EXPANSION  OF  INTERMEDIATE 
COEFFICIENTS,  The  SCALAR  VERSION  REQUIRES  19  OPERATIONS  PER  POINT 
while  this  vector  version  REQUIRES  23.  THESE  CAN  BE  SPLIT  ACROSS 
BOTH  PIPES,  however,  AND  ARE  PERFORMED  1N‘ VECTOR  MODE. 


ARGUMENTS!  (R  ARR(... 
NC  INTEGER 

A R ARR(NC) 

B R ARR(NC) 

C R ARR(NC) 

0 R ARR(NC) 

X R ARRtNC) 

SCA  R ARR(NC) 

SCB  R ARR(NC) 


) IS  USED  TO  denote  REAL  ARRAY  & DIMENSIONS) 


• EONS  IN  THE  TRIDIAGONAL  SYSTEM  I 
COEFFS  OF  X(I-l)  IN  I-TH  EQN.  I 
COEFFS  OF  X(I)  IN  I-TH  EQN.  I 
COEFFS  OF  X(Itl)  IN  I-TH  EQN.  I 
0(1)  X INHOMOGENEOUS  PART  OF  I-TH  EQN  1 
SOLUTION  VECTOR  SOUGHT  0 
SCRATCH  ARRAY  FROM  CALLING  PROGRAM  S 

SCRATCH  array  from  calling  program  S 


LANGUAGE  ANO  LIMITATIONS!  THIS  ROUTINE  IS  WRITTEN  IN  OOUBLE  PRE- 
CISION fortran  for  the  ASC  (64  BITS).  DIAGONAL  DOMINANCE  IS 
STRONGLY  RECOMMENOEO  FOR  GOOD  ERROR  CONTROL.  THIS  ROUTINE  IS  REC- 
OMMENDED WHEN  LONG  (NC  > IS)  TRIDIAGONAL  SYSTEMS  ARE  BEING  SOLVED. 
WHEN  CALLED  WITH  SHORTER  SYSTEMS  TRIOOV  CALLS  TRIDDS.  THE  SHORT 
RESIDUAL  SYSTEMS  GENERATED  BY  THE  VECTOR  FOLDING  OPERATIONS  ARE 
SOLVED  BY  CALLING  TRIOOS.  THE  MULTIPLE  TRIDIAGONAL  SOLVER  TRIOOM 
IS  BASED  ON  AN  EXPANSION  OF  TRIOOV. 

NOTE!  IN  TRIOOF  THERE  ARE  TWO  PARAMETERS  WHICH  CONTROL  THE  VEC- 
TOR SOLUTION.  LMAX  IS  THE  MAXIMUM  NUMBER  OF  FOLDING  OPERATIONS, 
LMAX  X u CURRENTLY.  IF  SYSTEMS  LONGER  THAN  ABOUT  2«*14  ARE  TO  BE 
SOLVED,  LMAX  MUST  BE  INCREASED.  WHEN  LMAX  IS  INCREASED,  LMIN  MUST 
BE  INCREASED  AT  LEAST  AS  MUCH. 


TRIOOv  has  no  ADDITIONAL  USER  ENTRIES  AND  CALLS  THE  SUBROUTINE 
TRIOOF  TO  perform  THE  VECTORIZED  FOLDING, 


C 

C 

C 


implicit  REAL*8  (A-M,0-Z) 

DIMENSION  A(NC),  6(NC),  C(NC),  0(NC),  X(NC),  SCA(NC),  SCB(NC) 


TRIOOF  manages  both  a ANO  ANEW,  ETC,  IN  THE  SAME 
VECTOR  HAZARDS. 

II  X NC/2 

CALL  TRIOOF  (NC,  A,  S,  C,  0,  X,  SCA,  SCS(l), 
1 SCB(2*II^1),  SCB(3*Iltl),  A,  B,  C,  0,  X, 

RETURN 


ARRAY  WITHOUT 


SCB(llYl), 

2*NC) 
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SUBROUTINE  TRIDOM  (NC«  A,  B,  C>  0,  X,  SCA,  SCU>  NS>  NO) 

CO 

CO  *•*••***•**•**•*•*•***•****•*••*• 

CO 

CO  TRIOON  (NC>  A>  B«  C«  0>  X,  SCA,  SCB>  NS>  NO)  CLASS. 

CO  ORIGINATOR  - J.  BORIS  COOE  7706,  NRL  JULY  1976 

CO 

CO  OESCRIPTION:  this  routine  solves  in  OOUBLt  PRECISION  THE  SET  OF 
CO  NS  parallel  TRIOIAGONAL  matrix  equations  max  s 0 hHICH  are  all  of 
CO  THE  SAME  LENGTH.  HERE  A(IS,I),  B(IS,I)  ANO  C(IS,I)  ARE  THE  ONLY 

CO  THREE  NON  ZERO  ELEMENTS  OF  ROW  I IN  SET  IS.  SPECIFICALLY  WE  SEEK 

CO  X(IS,I)  IN  THE  EQUATION 

CO 

CO  a(is,i)ax(is,i*i)  * b(is,i)ax(is,i)  * C(is,i)Ax(is,itn  * ons,n 

CO 

CO  FOR  is  s 1,  ...,  NS  ANO  I : 1,  ...,  NC.  FOR  PERIOOIC  BOUNDARY  CON- 

CO  OITIONS  X(IS,n  : X(IS,NCA1)  AND  X(IS,0)  : X(IS,NC).  THE  APERIODIC 

CO  CASE  FOLLOWS  FROM  SETTING  A(IS,1)  ANO  C(IS,NC)  TO  ZERO.  THE  ALGO- 

CD  RITHM  USED  HERE  FOLLOWS  FROM  REPEATED  CtCLIC  REDUCTION  (NRL  MEMO- 

CD  RANOUM  REPORT  3l4«,  OCT.  76).  THE  AUXILIARY  ROUTINE  TRIOOG  IS  USED 

CO  TO  AID  WITH  storage  FOR  THE  NX  COMPLIER.  ALL  ARRAYS  MUST  HAVE  AT 

CO  LEAST  2*NC  LOCATIONS  TO  ALLOW  FOR  EXPANSION  OF  INTERMEDIATE 
CO  COEFFICIENTS.  The  SCALAR  VERSION  REQUIRES  19  OPERATIONS  PER  POINT 

CO  while  this  vector  VERSION  REQUIRES  23.  THESE  CAN  BE  SPLIT  ACROSS 

CO  BOTH  PIPES,  however,  AND  ARE  PERFORMED  IN' VECTOR  MODE. 

CO 

CO  ARGuMENTSI  (R  ARR(,..)  IS  USED  TO  DENOTE  REAL  ARRAY  ( DIMENSIONS) 
CO  NC  INTEGER  » EONS  IN  THE  TRIOIAGONAL  SYSTEM  I 

CD  A R ARR(N0,2*NC)  COEFFS  OF  X(IS,I-1)  IN  I-TH  EQN,  I 

CO  B R ARRtND,2*NC)  COEFFS  OF  X(IS,I)  IN  I-TH  EQN.  I 

CD  C R ARR(NP,2*NC)  COEFFS  OF  X(1S,I+1)  IN  I-TH  EQN,  I 

CO  0 R ARR(N0,2*NC)  0(IS,I)  s INHOMOGENEOUS  PART  OF  EQN  I I 

•CO  X R ARRtN0,2*NC)  SOLUTION  VECTOR  SOUGHT  X(IS,I)  0 

CO  SCA  R ARR(N0,2*NC)  SCRATCH  ARRAY  FROM  CALLING  PROGRAM  S 

CD  SCB  R ARRtN0,2*NC)  SCRATCH  ARRAY  FROM  CALLING  PROGRAM  S 

CO  NS  integer  • SETS  OF  TRIOIAGONAL  SYSTEMS  I 

CO  NO  integer  DIMENSION  OF  FIRST  INDEX  (,GE.  NS)  I 

CO 

CO  LANGUAGE  ANO  LlMIjATlONSI  THIS  ROUTINE  IS  WRITTEN  IN  OOuBLE  PRE- 

CO  CISION  FORTRAN  FOR  THE  ASC  (64  BITS).  DIAGONAL  DOMINANCE  IS 

CD  STRONGLY  RECOMMENDED  FOR  GOOD  ERROR  CONTROL.  THIS  ROUTINE  IS 
CO  PARTICULARLY  EFFICIENT  WHEN  MANY  LONG  TRIOIAGONAL  SYSTEMS  (NS>12) 

CO  HAVE  TO  BE  SOLVED  IN  PARALLEL.  TRIOOR  WOULD  BE  BETTER  FOR  MANY 

CO  SHORT  TRIOIAGONAL  SYSTEMS.  TRIOOH  CALLS  TRIOOR  TO  SOLVE  THE  SERIES 

CO  OF  SHORT  RESIDUAL  SYSTEMS  BEING  GENERATED  BY  THE  VECTOR  FOLDING 

CO  OPERATIONS.  THE  FX  COMPILER  SHOULD  BE  USED  TO  AVOID  A COMPILER  BUG 
CD  NOTE!  IN  TRIOOG  THERE  ARE  TWO  PARAMETERS  WMICH  CONTROL  THE  VEC- 
CO  TOR  solution,  LMAX  IS  THE  MAXIMUM  NUMBER  OF  FOLDING  OPERATIONS. 

CO  LMAX  a 10  CURRENTLY.  IF  SYSTEMS  LONGER  THAN  ABOUT  2**13  ARE  TO  BE 

CD  SOLVED,  LMAX  MUST  BE  INCREASED.  WHEN  LMAX  IS  INCREASED,  LMIN  MUST 
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o r>  o 


CO  BE  INCREASED  AT  LEAST  AS  MUCH. 

CO 

CO  TRIOOM  HAS  NO  ADDITIONAL  USER  ENTRIES  AND  CALLS  THE  SUBROUTINE 
CO  TRIOOG  TO  PERFORM  THE  VECTORIZED  FOLDING, 

CO 

CO  ••*••*•*••***•*•*•••«*«*•**•••••• 

c 

implicit  REAL»8  (A-H,0-Z) 

DIMENSION  ACNO^NC)/  B(ND,NC}«  C(ND,NC)/  D(NO>NC)>  X(NOfNC) 
DIMENSION  SCACNOfNC),  SCbCND.NC) 

TRIOOG  manages  BOTH  A AND  ANEW,  ETC.  IN  THE  SAME  ARRAY  WITHOUT 
VECTOR  HAZARDS. 

II  s NC/2 

CALL  TRIOOG  (NC.  A.  B.  C.  0,  X.  SCA.  SCBU.D.  SCB(1.II«1). 

1 SCB(lf2«II«n.  SCBa.3*IItl).  A.  B,  C.  0.  X.  2*NC.  NS.  NO) 

RETURN 

End 
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subroutine  TRIOOR  lUCt  i,  8«  Cl  Of  X,  R6f  AiNVQf  NSf  NO) 

CO 

CO  •*.**•.»*****.********«.*.*,*,,,* 

CD 

CO  TRIOOR  (NCf  Af  Sf  Cf  Of  X,  RB,  AINVQf  NSf  NO)  CLASS. 

CO  ORIGINATOR  • J.  BORIS  CODE  TTObf  NrL  JULY  1976 

CD 

CO  OESCRIPTIONj  THIS  ROUTINE  SOLVES  IN  DOUBLE  PRECISION  THE  SET  OF 
CO  NS  parallel  TRIOIAGONAL  MATRIX  equations  H*X  : 0 WHICH  ARE  ALL  OF 
CO  THE  SAME  LENGTH.  HERE  A(lSfI)f  B(ISfl)  AND  C(ISfl)  ARC  THE  ONLY 

CO  THREE  NON  ZERO  ELEMENTS  OF  ROW  I IN  SET  IS.  SPECIFICALLY  WE  SEEK 

CD  X(ISfl)  IN  THE  EQUATION 

CD 

CO  A(is,i)*x(iSf i-i)  * B(iSf i)*x(iSf n * cdSf i)*xciSf ui)  > ociSfi) 
CO 

CO  FOR  IS  S If  ...f  NS  AND  I : If  ...f  NC . FOR  PERIODIC  BOUNDARY  CON- 
GO DITIONS  X(ISfl)  3 X(ISfNCtl)  AND  X(ISfO)  s X(ISfNC).  THE  APERIODIC 

CD  CASE  FOLLOWS  FROM  SETTING  A(ISfl)  ANO  C(ISfNC)  TO  ZERO.  THE  SCALAR 

CO  ALGORITHM  EMPLOYED  HERE  FOLLOWS  FROM  SEPARATING  M INTO  A PURELY 

CO  TRIOIAGONAL  PART  ANO  A VECTOR  Q IN  COLUMN  NC  WITH  ONE  AT  ROW  NC . 

CO  SINCE  NS  CASES  ARE  BEING  SOLVED  IN  PARALLEL.  VECTORIZaTION  ON  THE 

CO  CASES  IS  PROVIDED.  NOTE  THAT  ALL  ARRAYS  OF  INPUT  COEFFICIENTS  ANO 

CO  DATA  MAY  have  TO  BE  TRANSPOSED  TO  ENSURE  THAT  THE  CORRESPONDING 
CO  elements  of  data  in  the  ns  SYSTEMS  ARE  CONTIGUOUS.  NO  IS  THE  INC- 

CO  REMENT  IN  MEMORY  TO  THE  NEXT  ELEMENT  OF'THE  SAME  SYSTEM. 

CO  THIS  ROUTINE  IS  PARTICULARLY  EFFICIENT  WHEN  MANY  TRIOIAGONAL 

CO  SYSTEMS  (NS  > 20)  HAVE  TO  BE  SOLVED  IN  PARALLEL.  TRIOOS  MOULD  BE 

CO  better  for  only  a few  short  systems,  the  primary  purpose  of  TRIOOR 

CO  IS  TO  complete  the  folding  calculation  IN  trioom, 

CO 

CO  ARGUMENTS!  (R  ARR(...)  IS  USED  TO  DENOTE  REAL  ARRAY  & DIMENSIONS) 


CD 

NC 

INTEGER 

» EONS  IN 

THE  TRIOIAGONAL  SYSTEM 

I 

CO 

A 

R 

ARR(NOfNC) 

COEFFS  OF 

X(ISfl-l)  IN  I-TH  EQN. 

I 

CO 

B 

R 

ARR(NDfNC) 

COEFFS  OF 

X(ISfl)  IN  I-TM  EQN, 

I 

CD 

C 

R 

ARR(NOf NC) 

COEFFS  OF 

X(IS,I+l)  IN  I-TH  EQN, 

I 

CO 

D 

R 

ARR(NOfNC) 

D(ISfl)  = 

INHOMOGENEOUS  PART  OF  EQ  I 

I 

CO 

X 

R 

ARR(NOfNC) 

SOLUTION 

VECTOR  sought 

0 

CO 

RB 

R 

ARR(NOfNC) 

SCRATCH  array  from  CALLING  PROGRAM 

s 

CO 

AINVO 

R 

ARRtNDfNC) 

SCRATCH  ARRAY  FROM  CALLING  PROGRAM 

s 

CD 

NS 

INTEGER 

• SETS  OF 

TRIOIAGONAL  SYSTEMS 

I 

CD 

CD 

NO 

INTEGER 

DIMENSION 

OF  FIRST  INDEX  (,GE.  NS) 

1 

CO  language  AND  LIMITATIONS!  THIS  ROUTINE  IS  WRITTEN  IN  DOUBLE  PRE' 

CO  CISION  FORTRAN  FOR  THE  ASC  (64  BITS).  DIAGONAL  DOMINANCE  IS 

CO  STRONGLY  RECOMHENOEO  FOR  GOOD  ERROR  CONTROL.  AS  WRITTEN,  NS  < 202 
CD  IS  REQUIRED  BECAUSE  OF  THE  PARAMETER  NPT  BELOW. 

CD 

CO  TRIOOR  MAS  NO  ADDITIONAL  USER  ENTRIES  AND  CALLS  NO  AUXILIARY  SuB* 
CO  ROUTINES. 

CO 
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CD 

C 


C 


C 


C 

C 


C 


C 

C 


c 


c 


implicit  REAL*0  (A-H.O-Z) 
parameter  NPT  « 202 

DIMENSION  A(ND>NC)>  a(NO»NC)>  C(ND>NC}>  D(NOfNC),  X(ND>NC) 
DIMENSION  R8(N0>NC)>  AINVO(ND>NC) f RCAQ(NPT) 

NM  a NC  • I 

DO  10  I s If  NC 
DO  10  IS  = If  NS 
IP  (B(ISfl)  .NE.  0.0)  GO  TO  10 
B(lSfl)'=  l.OE-70 
10  CONTINUE 

DO  12  IS  « If  NS 
12  RBIISfl)  a 1.000/B(IS.  1) 

DO  IS  I a 2f  NH 
DO  IS  IS  a If  NS 

IS  RB(lS,I)al.D+0/  (B(ISfl)  • A(ISfl)  * C(ISfl-l)  * RB(ISfl-l)) 

INVERT  TRIOIA60NAL  PART  ON  Q 
DO  17  IS  a If  NS 

17  AINVQ(ISfl)  a A(ISfl)  * RB(ISfl) 

CO  20  I a 2,  NM' 

DO  20  IS  a If  NS 

20  AINVQ(ISfl)  a . A(lSfl)  • AINVQ (ISf !•! } * RB(ISfl) 

DO  2S  IS  a 1,  NS 

2S  AINVQ(lSfNM)  a AINVQ(ISfNM)  * CCISfNM)  • RB(ISfNM) 

DO  30  J a 2f  NM 
I a NC  - J 
DO  30' IS  a 1,  NS 

30  AINVO(ISfl)  a AINVOdSfl)  - C(ISfl)  « RB(ISfl)  • AlNVQdSfl^l) 
DO  3SIS  a If  NS  ■ 

3S  RCAQdS)  a 1 . ODO/ CB dSf NC } - C(ISfNC)  • AINVO(ISfl) 

1 - AdSfNC)  * AlNVQ(ISfNM)} 

END  precomputation.  the  following  are  all  the  CALCS  INVOLVING  0 
DO  37  IS  a If  NS 

37  XdSfl)  a O(ISfl)  a RBdSfl) 

DO  40  I a 2i  NM 
DO  40  IS  a If  NS 

40  XdSfl)  a (DdSfl)  - AdSfl)  • XdSfI-1))  * RB(ISfl) 

DO  SO  J a 2f  NM 
I a NC  - J 
00  SO  IS  a If  NS 

SO  XdSfl)  a X(ISfl)  - C(ISfl)  • RBdSfl)  * X(ISflTl) 

00  SS  IS  a If  NS 


5S  X(ISfNC)  a (O(ISfNC)  - C(ISfNC)  * X(ISfl)  - A (I S f NC ) *X ( I S f NM ) ) 
1 * RCAQ(IS) 

DO  bO  I a If  nM 

00  60  IS  a If  NS 

60  XdSfl)  a XdSfl)  • XdSfNC)  * AINVOdSfl) 

RETURN 

ENO 


18 


CO 

CO 

CO 

CO 

CO 

CO 

CO 

CD 

CD 

CO 

CO 

c 

c 


SUBROUTINE  TRIDDF  (NC,  A,  B,  C«  0,  X,  RB,  SCA,  SCB«  SCC>  SCO, 
1 ANEW,  BNEM,  CNEW,  ONEM,  XNEM,  LENGTH) 


TRIDOF  ( — — .— ) CLASS. 

ORIGINATOR  - J,  BORIS  CODE  7706,  NRL  JULY  1976 

TRIDOF  IS  AN  internally  USED  AUXILIARY  ROUTINE  TO  TRIDOV  AND  NEED 
NOT  BE  REFERENCED  DIRECTLY  BY  THE  USER,  TRIDOF  IS  SEPARATE  FROM 
TRIDOV  SO  THAT  A & ANEW,  B & BNEW,  ETC.  CAN  BE  EQUI V ALENCED  ' V I A 
The  ARGUMENT  LISTS. 


NOTE  THAT  LMIN  SHOULD  EXCEED  LMAX, 
implicit  REAL*B  (A-M,0-Z) 

parameter  lmax  » n,  lmin  s i4,  lmaxp  = lmax  ♦ i 

DIMENSION  NL(LMAXP),  ML(LMAXP),  lO(LMAXP),  JOCLMAXP) 

DIMENSION  A(LENGTH),B(LENGTH),C(LENGTH),  D(LENGTM),  X(LENGTM) 
DIMENSION  RB(LENGTH),  ANEW (LENGTH) , BNEW  (LENGTH ) ,CNEW (LENGTH ) 
DIMENSION  ONEW(LENGTH) , XNEW(LENGTh) 

DIMENSION  SCA(NC),  SCB(NC),  SCC(NC),  SCD(NC) 

FOLD  the  BEJESUS  OUT  OF  THE  COEFFICIENTS. 

N = NC 


LEVEL  s 0 
I0(LEVEL»1)  = N 
JO(LEVEL^l)  < 0 
DO  500  Lai,  LMAX 
II  = IO(L) 

J1  a JO(Ll 
NL(L)  a’N 

IF  (N  ,LE.  LMIN)  GO  TO  600 

M a (Ntl)/2 

HH  X M • 1 

ML(L)  a M 

DO  200  II  a 1,  M 

I a II  ♦ II 

J a 2.11  t J1  - 1 

BNEW(I)  a B(J) 

ONEM(I)  a D(J) 

200  RB(I)  a 1.0/B(J^1) 

DO  210  II  a 1,  MM 
I a II  . II 
J a 2*11  * J1 
SCA(Il)  a A(J*1)*RB(I) 
SCC(II)  a C(J-1)*RB(I) 
ANEW(Itl)  a A(J)*3CA(1I) 
CNEW(I)  a C(J)*SCC(in 
ANEW(I.l)  a -ANEW(I.l) 
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CNEW(I)  = -CNEH(I) 

SCSdl)  3 C(J}*SCA(II) 

SCOdl)  3 D(J)*SCA(d} 

BNEWd  + 1)  3 BNEWd  + 1)  - SCB(d) 

ONEWCI  + l)  = ONEWd  + l)  - SCDdI) 

SCBdl)  3 A(J)*SCC(n) 

SCO(d)  3 OtJ)*SCC(d) 

BNEM(I)  3 BNEWd)  • SCBCII) 

210  ONEHd)  3 DNEW(l)  - SCDCin 

CAUL  BOTCHO  ANEW,  BNEW,  CNEW,  DNEw,  XNEW) 

IF  (2*M  ,NE.  N)  GO  TO  220 
TA  3 A(JUn/B(JltN) 

TC  3 C(JUN-l)/B(Jl+N) 

ANEWdltl)  3 . TA*A(J1  + N) 

BNEWdUd  3 BNEWdltl)  • TA*C(JltN) 

ONEWdltl)  3 DNEWdl  + 1)  - TA*OCJl+N) 

BNEWdl+Mj  3 BNEWdl  + M)  - TC*A(JltNj 
CNEWdl+Mj  3 -TC*C(JJtN) 

ONEWdl+M)  3 ONEWCn+H)  - TC*0(J1>NJ 
GO  TO  230 

220  ANEWdl^l)  3 ACJUl] 

CNEWdltMj  3 C(JUN) 

230  continue  ' 

JO(Ltl)  3 JO(L)  ♦ N 

IO(L  + n 3 IO(L)  t M 

N 3 M 

500  LEVEL  3 LEVEL  ♦ I 

C 

C USE  THE  SCALAR  VERSION  FOR  SHORTIES. 

600  CALL  TRIDDS  (N,  A(Jltl),  C(J1>1),  D(JUI), 

1 XCJltUf  SCA,’SCB) 

IF  (LEVEL  ,EQ.  0)  RETURN 
C 

C unfold  and  the  answer  appears  miraculously. 

00  900  LL  3 1,  level 
L 3 LEVEL  ♦ 1 - LL 
M 3 ML(L) 

N 3 NL(L) 

n 3 io(u 

J1  3 JO(L) 

DO  BOO  II  3 1,  M 

1 3 II  ^ 11 
J 3 2*11  ♦ Jl  - I 
X(J)  3 XNEW(I) 

MM  3 M « 1 

DO  810  II  3 1,  MM 
I 3 II  * II 
J 3 2*11  * Jt 
SCAdI)  3 A(J}3XNEW(J*1) 

SCCCin  3 C(J)*XNEW(J*1) 

SCBCII)  3 D(J)  - SCAdI) 

SCOCII)  3 SCBdl)  • SCC(II) 

X(J)  3 SCO(II)*RB(I) 

IF  (2*M  ,EQ.  N) 

X(Jt*N)  3 (0(J1*N)  - A(JltN)*X(JI*N-l)  - C(JltN)*X(Jltl)) 
/B(JltN) 

CONTINUE 
RETURN 

END 


800 


810 


900 
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SU8R8UTINE  TRIODG  (NC,  A,  0,  C.  0,  X,  Rbr  SCA,  SCB,  SCC,  SCO, 
1 ANEM,  8NEM,  CNEW,  OnEM,  XNEM,  length,  ns,  NO) 


CO 

CO  ********************************* 

CO 

CO  TRIOOG  ( — — — — ) CLASS. 


CO  ORIGINATOR  - J.  BORIS  COOE  7706,  NRL  JULT  1976 

CO  TRIOOG  IS  AN  INTERNALLY  USEO  AUXILIARY  ROUTINE  TO  TRIOOM  ANO  NEED 

CO  NOT  BE  REFERENCEO  OIRECTLY  BY  THE  USER.  TRIOOG  IS  SEPARATE  FROM 

CO  TRIOOM  SO  THAT  A & ANEW,  B & BNEW,  ETC.  CAN  BE  EOUIVALENCEO  VIA 

CO  THE  ARGUMENT  LISTS, 

CO 

CO  *.**»»***.*.*.**».»***»****».***. 

C 

C NOTE  THAT  LMIN  SHOULO  EXCEEO  LMAX. 

implicit  REAL*8  (A-H,0-Z) 

parameter  LMAX  = 10,  LMIN  s u,  LMAXP  s LMAX  ♦ 1,  NPT  = 202 
DIMENSION  NL(LMAXP),  MLCLMAXP),  IO(LMAXP),  JO(LMAXP) 
DIMENSION  A(NO, LENGTH] , B (NO , LENGTH ) , C(NO, LENGTH) 

DIMENSION  D(NO, LENGTH) , X (NO , LENGTH ) , RB (NO, LENGTH) 

DIMENSION  SCA(NO,NC),  SC8(N0,NC),  SCC(ND,NC),  SCD(NO,NC) 
DIMENSION  ANEW(NO, LENGTH)  , BNEW (NO , LENGTH ) , CNEW  (NO , LENGTH ) 
DIMENSION  DNEWtNO, LENGTH) , XNEW (NO , LENGTH  ) , TA(NPT),  TC(NPT) 
C 

C FOLD  THE  BEJESUS  OUT  OF  THE  COEFFICIENTS. 

N 3 NC 
LEVEL  = 0 
lOCLEVELtl)  = N 
JOCLEVEL+l)  s 0 
00  SOU  LSI,  LMAX 
II  s I0(L) 

J1  = JO(LJ 
NL(U  s N 

IF  (N  ,LE.  LMIN)  GO  TO  600 
M 3 (N  + D/2 
MM  3 H • 1 
ML(U  * M 

00  200  II  = 1,  M 

1 s II  ♦ n 

J 3 2«II  ♦ Jl  - I 

DO  200  IS  3 1,  NS 
BNEWdSfl)  3 B(IS,J) 

ONEW(IS,I)  3 0(IS,J) 

200  RB(:S,I)  3 1.0/B(IS,J>1) 

00  210  II  3 i,  MM 

1 3 II  ♦ II 

J 3 2. II  ♦ Jl 

00  210  IS  3 1,  NS 

SCA(I3,II)  3 A(IS,Jtl)*RB(IS,I) 

SCC(IS,II)  3 C(IS,J>1)*RB(IS,I) 
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210 


215 

220 


225 

230 


500 


ANEW(IS>1«1)  a *(IS>J)*SCA(IS,in 
CNCW(1S>I)  a C(IS,J)*SCC(IS>II) 

ANEwdS/I  + n » -ANEi^dSf  I + l) 

CNEWdSd)  a -CNEMdSfl) 

SCBCISid]  a CdS«J}*SCA(ISdI} 

SCOCISrU)  a 0 d S • J ) *SC  A (I  S ,1 1 ) 

BNEWdSd^l)  a BNEWdSd^l)  • SCBClSdl) 

ONEWdSd^ll  a ONEMdS.I^l)  - SCOdSdl) 

SCBdSdl)  a AdS,j)*SCC(IS.II) 

SCDdSdl]  a OdS,J)  ASCCClS.in 
BNEMdS»n  a BNEMClSd)  • SCBdSdl) 

DNEMCIS,!}  a ONEHdSd)  • SCOCISdU 

CAUL  BOTCHO  CN,H,ll>Jl>  ANEW , BNEM , ' CNEM « ONEW,  XNEH) 

IF  {2*M  ,NE.  N)  CO  TO  220 

00  215  IS  a l,  NS 

TACIS)  a A(IS,JUl)/BdS>JUN) 

TCCIS)  a CdS,JltN-l)/BdS,Jl*N) 

ANEW(IS»Ilfl)  a - TA(IS)*AdS,JlAN) 

BNEW(IS,I1  + 1)  - TA(IS)*CdS,Jl+N) 
0NEH(IS,I1+1)  - TACISjaOCIS, Jl+N) 
BNEWCIS,  Il+M)  - TCdSj*AdS,Jl+N) 
-TCCIS)*C(IS,J1+N) 

0NEW(1S,I1+M)  - TCdS)»OdS,Jl+N) 


BNEW(IS,Iltl) 

0NEW(IS,I1+1) 

BNEWCISdltM) 

CNEWCIS.IITM) 

ONEWCIS, IltM) 

GO  TO  230 
CON.'INUC 
00  225  IS  a i, 
ANEWClSdltl)  a AdS,Jl4'l} 
CNEMCIStlltM)  a CClSfJltN) 
CONTINUE 

J0(L+1)  a JOCL)  ♦ N 
lOCL^l)  a lOCU  ♦ M 
N a M 

LEVEL  a LEVEL  ♦ 1 


NS 


C 

C USE  THE  SCALAR  REPEAT  VERSION  FOR  SHORTIES. 

oOO  CALL  TRIOOR  (N,  A(l,Jlfl),  B(1,JU1),  C(1«J1>1),  OCdJltl]* 
1 X(l,Jltl),  SCA,  SCB,  NS,  NO) 

IF  (LEVEL  .EQ.  0)  RETURN 
C 

C UNFOLD  ANO  the  ANSWER  APPEARS  MIRACULOUSLY. 

00  900  LL  a 1,  level 
L a level  ♦ 1 • LL 
M a MLCL) 

N a NL(L) 

II  a lOCL) 

J1  a JQCLl 

00  BOO  II  a 1,  M 

1 a II  ♦ n 

J a 2*11  ♦ J1  - 1 
00  800  IS  a 1,  NS 
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800  X(IS,J)  S XNEM(lS<n 

MM  = M - I 

00  810  II  3 1«  MM 

1 = II  ♦ II 

J = 2*11  ♦ J1 
00  810  IS  s 1>  NS 
SCA(IS.II)  a A(IS,J)*XNEM(IS>J*1) 
SCC(IS,II)  a C(IS>J)*XNEM(ISrJtl) 
SCB(IS.II)  a 0(IS«J)  • SCA(IS,in 
SC0(IS,IU  a SCB(IS,II)  • SCCCISfin 
810  X(IS>J)  a SCD(IS,II)*HB(IS>n 

11^  (2*M  ,NE.  N)  GO  TO  <>00 
00  820  IS  a 1,  NS 
TA(IS)  a A(IS,J1*N)*X(IS,J1*N-1) 
TC(IS)  a CtlS, J1*N)*X(IS,J1*1) 
X(lS,JltNJ  a 0(IS,J1+N)  - TACIS) 
XtIS»Jl*N)  a X(IS,J1+N)  • TC(IS) 

820  X(IS>J1*N)  a X Cl S , J 1 tN } /B ( I S , J 1 *N > 

900  continue 

return 

ENO 


Appendix  B 

Test  Program  for  Double  Precision  Solvers 
Including  Both  Single  and  Double  Precision  Results 


Because  the  single  aund  double  precision  test  programs  submitted 
to  SPLi  are  essentially  identical,  only  the  double  precision  driver  is 
included.  The  results  from  both  single  and  double  precision  tests  are 
included  for  comparison,  however.  The  timing  routine  SECOND 
(from  / ASG  UTIL,  * * * ) is  used  to  time  the  various  test  executions. 


/ ASG  PP0UI8»  USERCAT/O77/B50/COOYJ1/PPOLIB»USE=SmR 
/ ASG  UTIL»USERCAT/CC0109aS/TMS/UTlLlTY/OaJECT»USE=SHR 
/ FTN  FTNOPT»(K,M,Y,V,0,U),INsSYS.FIN,FTNTIME=JOOOO 
C PROGRAM  TO  TEST  THE  THREE  TRIOIAGONAL  SOLVERS  SPL  J976 

IMPLICIT  REAL*8  (A-H,  0-Z) 

PARAMETER  NPT  = 100,  NX  s 50,  NS  = 50,  NO  = 50 

DIMENSION  X(NPT),  XHOLD(NPT),  A(NPT),  B(NPT).  C(NPT),  O(NPT) 

DIMENSION  SCA(NPT),  SCB(NPT) 

DIMENSION  AM(ND,NPT),  BM(N0,NPT),  CM(ND,NPT),  0M(ND,NPT) 


DIMENSION 

XM(ND,NPT)  , 

SCAM(ND, 

NPT), 

SCBM(ND,nPT) 

1000 

format 

(' 

1 

NEW 

TEST 

CASE  ' 

, 12, 

• 

OF 

PERIODIC 

TRIDOS',  /) 

toot 

FORMAT 

(5X, 

13,  lPb012. 

4) 

1002 

format 

(' 

1 

NEW 

TEST 

CASE  • 

, 12, 

t 

OF 

APERIODIC 

TRIODS',/) 

1003 

format 

(' 

1 

NEW 

TEST 

CASE  • 

, 12, 

1 

OF 

PERIODIC 

TRIODV,  /) 

1004 

FORMAT 

(' 

1 

NEW 

TEST 

CASE  • 

, 12, 

1 

OF 

APERIODIC 

TRIDDV',/) 

1005 

FORMAT 

(5x, 

' I 

Ad) 

Bd) 

' 

Cd) 

D (I ) • , 

1 

1 

XHOLD(I) 

DELTAd) 

', 

/; 

1006 

format 

( ' 

1 

NEW 

TEST 

CASE  ' 

, 12, 

4 

OF 

PERIODIC 

TRIDDM',  /) 

1007 

FORMAT 

(' 

1 

NEW 

TEST 

CASE  ' 

, 12, 

OF 

APERIODIC 

TRIODM* ,/) 

CALL  INDUMP 
NCASE  3 1 
NXM  3 NX  • 1 

TWOPI  3 4.0D0*OATANn.0O0) 

C 

c START  the  loop  Over  test  cases. 

DO  500  ICASE  3 1,  NCASE 

AS  3 0FL0AT(ICASE)*Th0PI/DFL0AT(NX) 

AE  3 DFL0AT{ICASE)/DFL0AT(NX) 

C 

C set  uP  THE  initial  CONDITIONS  FOR  SINGLE  AND  MULTIPLE  CASES. 

DO  50  I 3 1,  NX 

XHOLD(I)  3 0SIN(DFL0AT(I)*AS)*DEXP(0FL0AT(I)*AE) 

A(I)  8 0.SD0*(1.0D0  • DSIN(DFL0AT(I)*AS)) 
cm  3 OEXP(-0FL0AT(I)*AE) 

50  8(1)  3 4,0  - Ad)  - C(I) 

DO  55  I 8 2,  NXM 

55  D(I)  3 A(I)*XM0LD(I-1)  ♦ B(I)*XH8L0(I)  ♦ C ( I ) *XH0LD ( I ♦ I ) 
D(l)  «:  Ad)*XH0LD(NX)  t Bd)*XH0LD(i)  t Cd)*XH0LD(2) 

P(NX)  8 AiHX)*XM0LD(NXM)  t 8(NX)«XH0LD(NX)  t C (NX ) *XH0L0 ( I ) 


24 


> 


00  60  I s 1.  NX 
00  60  IS  s h NS 
AM(IS*n  : *(I) 

BM(IS>I)  3 B(I) 

CM(ISa)  s CCI) 

60  0M(IS»I)  s 0(1) 

C 

C CALL  THE  TRIOIAGONAL  SOLVERS  ANO  PRINT  THt  RESULTS. 

NRITE  (6«  1000)  ICASE 
CALL  SECONO  (1.  DTIME) 

CALL  TRIOOS  (NX,  A,  B,  C,  0,  X,  SCA,  SCB) 

CALL  SECONO  (0,  OTIME) 

PRINT  1986,  DTIME 
MRITE  (6,  1005) 

00  65  Is  1,  NX 
65  X(I)  s X(I)  - XHOLO(I) 

NRITE  (6,  1001)  (I,  A(I),  BCD,  C(l),  0(1),  XHOLO(l),  X(I), 
1 1 s 1,  NX) 

WRITE  (6,  lOQi)  ICASE 
CALL  SECONO  (1,  OTIME) 

CALL  TRIDOV  (NX,  A,  B,  C,  0,  X,  SCA,  SCB) 

CALL  SECONO  (0,  OTIME) 

PRINT  1984,  OTIME 
NRITE  (6,  1005) 

1984  FORMAT  (lOX, 'TOTAL  OTIME  s ',020.5,'  SECONDS'/) 

00  70  I s l,  NX 
70  X(I)  s X(I)  - XHOLO(I) 

WRITE  (6,  1001)  (I,  A(l),  B(l),  C(I),  0(1),  XHOLO(I),  X(I), 
1 1 = 1,  NX) 

NRITE  (6,  1006)  ICASE 
CALL  SECONO  (1,  OTIME) 

CALL  TRIDOM  (NX,  AH,  BM,  CM,  OM,  XM,  SCAM,  SCBM,  NS,  NO) 
CALL  SECONO  (0,  OTIME) 

PRINT  1984,  OTIME 
WRITE  (6,  1005) 

DO  75  I = 1,  NX 

75  XM(l,I)  > XM(1,I)  - XMOLO(I) 

WRITE  (6,  1001)  (I,  AM(1,I),  BM(1,I),  CM(1,I),  0M(1,I), 

1 XHOLO(I),  XM(1,1),  1 s 1,  NX) 

C 

c for  the  aperiodic  case. 

A(l)  S 0.0 
CtNX)  s 0.0 
00  80  I s 2,  NXM 

80  0(1)  s A(I)*XM0L0(I-1)  ♦ B(I)>*XM8L0(I)  ♦ C ( I ) *XHOLO  ( I ♦!  ) 

0(1)  s A(1)*XHOLO(NX)  f B(1)*XH0LD(1)  t C ( 1 ) *XHOLO (2 ) 

OtNX)  s A(NX)*XHOLO(NXM)  t B (NX  ) «XHOLO (NX ) ♦ C (NX ) sXHOLO ( 1 ) 
DO  85  1 s'l,  NX 
DO  85  IS  s 1,  NS 
AM(IS,1)  s *(I) 
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65 

C 

C 


90 


95 


too 


500 

C 


CH(IS#1)  s C(I) 

0H(1S>I)  > 0(1) 

CALL  THE  TRIOIAGONAL  SOLVERS  ANO  PRINT  THE  RESULTS. 

WRITE  (6f  1002)  ICASE 
WRITE  (6.  1005) 

CALL  TRIDDS  (NX.  A.  B.  C,  0,  t,  SCA.  SCB) 

DO  90  1 « 1.  NX 
X(l)  X X(I)  - XHOLO(I) 

WRITE  (6.  1001)  (I.  A(I),  B(l).  C(l),  0(1).  XHOLO(I).  X(I), 
I X 1.  NX) 

WRITE  (6.  1004)  ICASE 
WRITE  (6.  1005) 

CALL  TRIDOV  (NX.  A.  B.  C.  0.  X.  SCA.  SCB) 

00  95  I X 1,  NX 

X(I)  X X(I)  • XHOLO(I) 

WRITE  (6.  1001)  (1.  Ad).  Bd).  Cd).  Od).  XHOLOd).  Xd). 
I X 1.  NX) 

WRITE  (6.  lOOT)  ICASE 
CALL  SECOND  (1.  OTIHE) 

CALL  TRIDOM  (NX.  AH.  BM.  CM.  DM.  XM.  SCAM.  SCBM.  NS.  NO) 
CALL  SECOND  (0>  OTIHE) 

PRINT  198tt.  OTIHE 
WRITE  (6.  1005) 

00  100  I X 1,  NX 

XMd.n  X XM(t,I)  - XHOLOd) 

WRITE  (b.  lOOl)  (I.  AH(l.I).  BH(l.I).  CH(t.I).  OH(l.I). 

XHOLOd).  XM(l.I).  1 X 1.  NX)  ' 

CONTINUE  ' 

STOP 

END 
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NEW  TEST  CASE  1 8F  PERIODIC  TRIOOS 


TOTAL  DTIME  * 0.109630-02  SECONDS 


I 

Acn 

8(1) 

cm 

om 

XMOLO(I)  OELTA(I) 

1 

4.68600-01 

2.55120 

00 

9.80200-01 

2,91290 

•01 

6,40590-02  -2.77560- 

17 

2 

4.37330-01 

2,60190 

00 

9,60790-01 

5,58590 

•01 

1.30450-01  -2,77560- 

17 

3 

4.06310-01 

2,65190 

00 

9,41760-01 

8,34370 

•01 

1,98970-01  4.16330- 

17 

4 

3.75660-01 

2.79120 

00 

9.23120-01 

1.11770 

00 

2,69400-01  -4,16330- 

17 

S 

3.45490-01 

2,74970 

00 

9.04640-01 

1.40770 

00 

3,41520-01  -1,24900- 

16 

6 

3,15940-01 

2,7.9710 

00 

8.86920-01 

1.70330 

00 

4,15060-01  -2.77560- 

17 

7 

2.87110-01 

2.84350 

00 

8,69360-01 

2.00330 

00 

4,89760-01  -8,32670- 

17 

8 

2.59120-01 

2.88670 

00 

8.52140-01 

2.30670 

00 

5,65340-01  -9.71450- 

17 

9 

2.32090-01 

2.93260 

00 

8.35270-01 

2.61220 

00 

6.41500-01  -1.24900- 

16 

10 

2.06110-01 

2.97520 

00 

8,18730-01 

2,91850 

00 

7.17920-01  -6.93890- 

17 

11 

1.81290-01 

3,01620 

00 

8.02520-01 

3.22420 

00 

7.94280-01  -2.08170- 

16 

12 

1,57730-01 

3.05560 

00 

7.86630-01 

3.52810 

00 

8,70230-01  -1.38780- 

16 

13 

1,35520-01 

3.09340 

00 

7,71050-01 

3,82860 

00 

9,45420-01  -1.66530- 

16 

14 

1.14740*01 

3.12950 

00 

7,55780-01 

4,12430 

00 

1,01950  00  -4.44090- 

16 

IS 

9.54920-02 

3,16370 

00 

7,40820-01 

4,41370 

00 

1.09210  00  -2.22040- 

16 

16 

7.78360-02 

3.19600 

00 

7.26150-01 

4.69520 

00 

1.16270  00  -2.22040- 

16 

17 

6,18470-02 

3.22640 

00 

7,11770-01 

4.96720 

00 

1.23120  00  -2.22040- 

16 

18 

4,75860-02 

3,25470 

00 

6.97660-01 

5.22830 

00 

1.29690  00  -2.22040- 

16 

19 

3,51120-02 

3.28100 

00 

6.83860-01 

5,47670 

00 

1.35960  00  -2.22040- 

16 

20 

2.44720-02 

3.30520 

00 

6.70320-01 

5,71090 

00 

1.41880  00  -2.22040- 

16 

21 

1,57080-02 

3.32720 

00 

6,57050-01 

5.92930 

00 

1.47410  00  0,00000 

00 

22 

8.85640-03 

3,34710 

00 

6.44040-01 

6,13020 

00 

1.52520  00  -2.22040- 

16 

23 

3.94260-03 

3.36460 

00 

6.31280-01 

6,31220 

00 

1,57160  00  -2.22040- 

16 

24 

9.86640-04 

3.38020 

00 

6, 18780-01 

6,47370 

00 

1.61290  00  -2.22040- 

16 

25 

0,00000  00 

3.39350 

00 

6.06530-01 

6.61310 

00 

1.64870  00  -2.22040- 

16 

26 

9,86640-04 

3.40450 

00 

5.94520-01 

6,72690 

00 

1.67870  00  -2.22040- 

16 

27 

3.94260-03 

3.41330 

00 

5.82750-01 

6,81980 

00 

1,70250  00  O.OOOOO 

00 

28 

8.65640-03 

3.41990 

00 

5.71210-01 

6,88440 

00 

1.71970  00  -2.22040- 

16 

29 

1. 57080-02 

3.42440 

00 

5.59900-01 

6.92120 

00 

1,72990  00  -2.22040- 

16 

30 

2,44720-02 

3.42670 

00 

5.48810-01 

6.92920 

00 

1,73290  00  -2.22040- 

16 

31 

3,51120-02 

3,42690 

00 

5,37940-01 

6.90700 

00 

1,72840  00  0.00000 

00 

32 

4.75860-02 

3,42510 

00 

5.27290-01 

6.85370 

00 

1,71600  00  -4.44090- 

16 

33 

6.18470-02 

3,42130 

00 

5.16850-01 

6.76820 

00 

1,69550  00  -2.22040- 

16 

34 

7,76360-02 

3,41550 

00 

5.06620-01 

6.64970 

00 

1.66660  00  -2,22040- 

16 

35 

9,54920-02 

3.40790 

00 

4.96590-01 

6.49730 

00 

1.62920  00  -2.22040- 

16 

36 

1.14740-01 

3.39850 

00 

4.86750-01 

6.31040 

00 

1,58300  00  -2.22040- 

16 

37 

1.35520-01 

3.38740 

00 

4,77110-01 

6.08840 

00 

1.52790  00  -2,22040- 

16 

38 

1,57730-01 

3,37460 

00 

4.67670-01 

5.83090 

00 

1.46380  00  -2.22040- 

16 

39 

1,81290-01 

3.36030 

00 

4.58410-01 

5.53760 

00 

1.39050  00  -2.22040- 

16 

40 

2,06110-01 

3.34460 

00 

4.49330-01 

5,20840 

00 

1.30810  00  -2.22040- 

16 

41 

2.32090-01 

3.32750 

00 

4.40430-01 

4.84330 

00 

1.21660  00  -2.22040- 

16 

42 

2.59120-01 

3.30920 

00 

4,31710-01 

4.44240 

00 

1.11590  00  -2,22040- 

16 

43 

2,87110-01 

3,28970 

00 

4.23160-01 

4.00600 

00 

1.00620  00  -2.22040- 

16 

44 

3,15940-01 

3.26930 

00 

4.14780-01 

3.53470 

00 

8.87510-01  -5.55110- 

17 

45 

3.45490-01 

3.24790 

00 

4,06570-01 

3.02900 

00 

7,60060-01  -9,71450- 

17 

46 

3.75660-01 

3.22580 

00 

3.98520-01 

2.48970 

00 

6,24040-01  -8,32670- 

17 

47 

0.06310-01 

3,20310 

00 

3,90630-01 

1,91790 

00 

4.79690-01  -4.16330- 

17 

48 

4.37330-01 

3,17980 

00 

3.82890-01 

1.31470 

00 

3.27330-01  -1.66530- 

16 

49 

4,68600-01 

3.15610 

00 

3,75310-01 

6,81410 

• 01 

i. 67300-01  1.38780- 

17 

SO 

5.00000-01 

3.13210 

00 

3.67880-01 

1,07220 

•01 

1.42220-15  -1.52030- 

18 

27 


NEW  TEST  CASE  I OF  PERIODIC  TRIOOV 


TOTAL  OTIME  » 0.799680-03  SECONDS 


1 

A(l) 

8(1) 

C(l) 

0(1) 

XHOLO(l) 

OELTA(I) 

1 

4.68600-01 

2.55120 

00 

9.80200-01 

2.91290 

-01 

6.40590-02 

0.00000 

00 

2 

4.37330-01 

2.60190 

00 

9.60790-01 

5.58590-01 

1.30450-01 

-1.38780- 

17 

3 

4.06310-01 

2.65190 

00 

9,41760-01 

8.34370' 

•01 

1,98970-01 

1.38780- 

17 

4 

3.75660-01 

2,70120 

00 

9.23120-01 

1.11770 

00 

2.69400-01 

-8.32670- 

17 

5 

3.45490-01 

2,74970 

00 

9.04840-01 

1.40770 

00 

3.41520-01 

-0.32670- 

17 

6 

3.15940-01 

2.79710 

00 

8.86920-01 

1.70330 

00 

4.15060-01 

-9.71450- 

17 

7 

2.87110-01 

2.64350 

00 

8.69360-01 

2.00330 

00 

4.09760-01 

-2,22040- 

16 

8 

2.59120-01 

2.88870 

00 

8.52140-01 

2.30670 

00 

5.65340-01 

-1,52660- 

16 

9 

2.32090-01 

2.93260 

00 

8.35270-01 

2.61220 

00 

6.41500-01 

-2,35920- 

16 

10 

2.06110-01 

2.97520 

00 

8,18730-01 

2.91850 

00 

7.17920-01 

-6.93890- 

17 

11 

1.81290-01 

3.01620 

00 

8.02520-01 

3.22420 

00 

7.94280-01 

-2.63660- 

16 

12 

1.57730-01 

3.05560 

00 

7,86630-01 

3.52810 

00 

0.70230-01 

-1,24900- 

16 

13 

1.35520-01 

3.09340 

00 

7,71050-01 

3.82860 

00 

9.45420-01 

0.32670- 

17 

14 

1.14740-01 

3.12950 

00 

7,55780-01 

4.12430 

00 

1.01950  00 

-2,22040- 

16 

15 

9.54920-02 

3.16370 

00 

7.40820-01 

4,41370 

00 

1.09210  00 

-2.22040- 

16 

16 

7.78360-02 

3.19600 

00 

7.26150-01 

4.69520 

00 

1.16270  00 

-2.22040- 

16 

17 

6.18470-02 

3,22640 

00 

7.117'':-01 

4.96720 

00 

1.23120  00 

-4.44090- 

16 

18 

4.75860-02 

3.25470 

00 

6.97o60-01 

5.22830 

00 

1.29690  00 

-2.22040- 

16 

19 

3.51120-02 

3.28100 

00 

6.83060-01 

5.47670 

oc 

1.35960  00 

0,00000 

00 

20 

2.44720-02 

3.30520 

00 

6.70320-01 

5,71090 

00 

1.41880  00 

-2.22040- 

16 

21 

1.57080-02 

3.32720 

00 

6.57050-01 

5.92930 

00 

1.47410  00 

-2.22040- 

16 

22 

8,85640-03 

3.34710 

00 

6.44040-01 

6.13020 

00 

1.52520  00 

-2.22040- 

16 

23 

3.94260-03 

3.36480 

00 

6.31280-01 

6.3122D 

00 

1,57160  00 

-2.22040- 

16 

24 

0.86640-04 

3.38020 

00 

6.10700-01 

6.47370 

00 

1.61290  00 

-2.22040- 

16 

25 

0.00000  00 

3,39350 

00 

6.06530-01 

6.61310 

00 

1.64870  00 

-2,22040- 

16 

26 

9,86640-04 

3,40450 

00 

5.94520-01 

6.72890 

00 

1.67870  00 

-2.22040- 

16 

27 

3.94260-03 

3,41330 

00 

5,82750-01 

6,81980 

00 

1.70250  00 

-2.22040- 

16 

28 

8,85640-03 

3.41990 

00 

5.71210-01 

6,88440 

00 

1,71970  00 

-2.22040- 

16 

29 

1.57080-02 

3.42440 

00 

5,59900-01 

6.92120 

00 

1,72990  00 

-2.22040- 

16 

30 

2.44720-02 

3.42670 

00 

5,48010-01 

6.92920 

00 

1.73290  00 

-2.22040- 

16 

31 

3,51120-02 

3,42690 

00 

5,37940-01 

6.90700 

00 

1,72840  00 

-2.22040- 

16 

32 

4,75860-02 

3.42510 

00 

5,27290-01 

6.85370 

00 

1.71600  00 

-2,22040- 

16 

33 

6.18470-02 

3.42130 

00 

5.16050-01 

6.76820 

00 

1.69550  00 

-4,44090- 

16 

34 

7,78360-02 

3.41550 

00 

5.06620-01 

6.64970 

00 

1,66660  00 

-2.22040- 

16 

35 

9,54920-02 

3.40790 

00 

4.96590-01 

6,49730 

00 

1.62920  00 

-2.22040- 

16 

36 

1.14740-01 

3.39850 

00 

4.86750-01 

6.31040 

00 

1.58300  00 

-2.22040- 

16 

37 

1.35520-01 

3.38740 

00 

4.77110-01 

6.0884O 

00 

1.52790  00 

-4.44090- 

16 

38 

1.57730-01 

3.37460 

00 

4,67670-01 

5.83090 

00 

1.46380  00 

-2.22040- 

16 

39 

1,81290-01 

3.36030 

00 

4.58410-01 

5.53760 

00 

1.39050  00 

-2.22040- 

16 

00 

2.06110-01 

3.34460 

00 

4.49330-01 

5.20840 

00 

1.30810  00 

-2.22040- 

16 

01 

2.32090-01 

3.32750 

00 

4.40430-01 

4.84330 

00 

1.21660  00 

-2.22040- 

16 

02 

2.59120-01 

3.30920 

00 

4.31710-01 

4.44240 

00 

1.11590  00 

-2,22040- 

16 

43 

2.87110-01 

3,28970 

00 

4,23160-01 

4.00600 

00 

1.00620  00 

-2.22040- 

16 

44 

3.15940-01 

3.26930 

00 

4,14700-01 

3,53470 

00 

8.87510-01 

-1 .52660- 

16 

05 

3.45490-01 

3.24790 

00 

4.06570-01 

3.02900 

00 

7.60060-01 

-2,08170- 

16 

46 

3.75660-01 

3.22580 

00 

3.98520-01 

2.48970 

00 

6,24040-01 

-9,71450- 

17 

47 

4.06310-01 

3.20310 

00 

3.90630-01 

1.91790 

00 

4,79690-01 

• 1 ,80410- 

16 

48 

4.37330-01 

3.17980 

00 

3.82890-01 

1.31470 

00 

3.27330-01 

-1,94290- 

16 

09 

4.68600-01 

3.15610 

00 

3.75310-01 

6.81410 

•01 

i, 67300-01 

1,38780 >17 

50 

5.00000-01 

3.13210 

00 

3.67880-01 

1,07220 

-01 

1.42220-15 

1,34140- 

19 

28 


NEW  TEST  CASE  I 8F  PERIODIC  TRIDOM 


TOTAL  OTIME  s 0.154630'01  SECONDS 


I 

Ad) 

am 

cm 

0(1} 

XHOLOd} 

OELTAd) 

1 

4.6060D-01 

2.55120 

00 

9,80200-01 

2,91290. 

01 

6.40590-02 

0,00000 

00 

2 

4.3733D-01 

2.60190 

00 

9,60790-01 

5.58590- 

01 

1,30450-01 

-1.38780 

17 

3 

4.0631D-01 

2.65190 

00 

9,41760-01 

8,34370- 

01 

1.98970-01 

1.38780 

17 

4 

3.75660-01 

2.70120 

00 

9.23120-01 

1.11770 

00 

2.69400-01 

-6.93890 

17 

S 

3.45490-01 

2.74970 

00 

9,04840-01 

1,40770 

00 

3.41520-01 

-9.71450 

17 

6 

3.15940-01 

2.79710 

00 

8.86920-01 

1.70330 

00 

4.15060-01 

-9.71450 

17 

7 

2.87110-01 

2,84350 

00 

8.69360-01 

2.00330 

00 

4.89760-01 

-2.22040 

16 

S 

2.59120-01 

2,88870 

00 

8.52140-01 

2.30670 

00 

5.65340-01 

-1.52660 

16 

4 

2.32090-01 

2.93260 

00 

8.35270-01 

2,61220 

00 

6.41500-01 

-2.77560 

16 

10 

2. 06110-01 

2,97520 

00 

8,i8730-01 

2.91850 

00 

7.17920-01 

-6.93890 

17 

11 

1.81290-01 

3.01620 

00 

8.02520-01 

3.22420 

00 

7.94280-01 

-2.63680 

16 

12 

1.57730-01 

3.05560 

00 

7.86630-01 

3,52810 

00 

8.70230-01 

-1.24900 

16 

13 

1.35520-01 

3.09340 

00 

7.71050-01 

3.82860 

00 

9,45420-01 

6.93890 

17 

14 

1.14740-01 

3.12950 

00 

7,55780-01 

4.12430 

00 

1,01950  00 

-2.22040 

16 

IS 

9.54920-02 

3,16370 

00 

7.40820-01 

4,41370 

00 

1.09210  00 

-2.22040 

16 

16 

7.78360-02 

3.19600 

00 

7.26150-01 

4,69520 

00 

1.16270  00 

-2.22040 

16 

17 

6.18470-02 

3,22640 

00 

7,11770-01 

4,96720 

00 

1.23120  00 

-2.22040 

16 

IS 

4.75860-02 

3,25470 

00 

6.97680-01 

5.22830 

00 

1,29690  00 

-2.22040 

16 

19 

3.51120-02 

3,28100 

00 

6.83860-01 

5,47670 

00 

1.35960  00 

0.00000 

00 

20 

2.44720-02 

3.30520 

00 

6,70320-01 

5.7109O 

00 

1.41880  00 

-2.22040 

16 

21 

- 1.57080-02 

3,32720 

00 

6,57050-01 

5.92930 

00 

1.47410  00 

-2.22040 

16 

22 

8.85640-03 

3.34710 

00 

6.44040-01 

6.13020 

00 

1.52520  00 

-2.22040 

16 

23 

3.94260-03 

3.36480 

00 

6.31280-01 

6.31220 

00 

1.57160  00 

-2.22040 

16 

24 

9.86640-04 

3.38020 

00 

6.18780-01 

6.47370 

00 

1.61290  00 

-2.22040 

16 

2S 

0.00000  00 

3,39350 

00 

6.06530-01 

6.61310 

00 

1.64870  00 

0,00000 

00 

26 

9.86640-04 

3.40450 

00 

5.94520-01 

6.72890 

00 

1,67870  00 

-2.22040 

16 

27 

3.94260-03 

3.41330 

00 

5.82750-01 

6.81980 

00 

1.70250  00 

•2.22040 

16 

2S 

8.85640-03 

3,41990 

00 

5.71210-01 

6.88440 

00 

1.71970  00 

-2.22040 

16 

29 

1.57080-02 

3.42440 

00 

5.59900-01 

6.92120 

00 

1.72990  00 

0.00000 

00 

30 

2.44720-02 

3.42670 

00 

5.48810-01 

6,92920 

00 

1,73290  00 

-2.22040 

16 

31 

3.51120-02 

3.42690 

00 

5.37940-01 

6.90700 

00 

1.72840  00 

-2.22040 

16 

32 

4.75860-02 

3.42510 

00 

5,27290-01 

6,85370 

00 

1.71600  00 

-2.22040 

16 

33 

6.18470-02 

3.42130 

00 

5,16850-01 

6.76820 

00 

1,69550  00 

-4.44090 

16 

34 

7.78360-02 

3.41550 

00 

5.0662D-01 

6.64970 

00 

1,66660  00 

-2.22040 

16 

3S 

9.54920-02 

3,40790 

00 

4.96590-01 

6,49730 

00 

1.62920  00 

-2.22040 

16 

36 

1.14740-01 

3.39850 

00 

4.86750-01 

6.31040 

00 

1.58300  00 

-2,22040 

16 

37 

1.35520-01 

3.38740 

00 

4.77110-01 

6,00040 

00 

1.52790  00 

-2.22040 

16 

38 

1.57730-01 

3.37460 

00 

4,67670-01 

5,83090 

00 

1.46380  00 

-2.22040 

16 

39 

1.81290-01 

3.36030 

00 

4.58410-01 

5.53760 

00 

1.39050  00 

-2.22040 

16 

do 

2,06110-01 

3.34460 

00 

4.49330-01 

5.20840 

00 

1.30810  00 

-2,22040 

16 

41 

2.32090-01 

3.32750 

00 

4,40430-01 

4.84330 

00 

1.21660  00 

-2.22040 

16 

42 

2.59120-01 

3.30920 

00 

4.31710-01 

4,44240 

00 

1.11590  00 

-2.22040 

16 

43 

2.87110-01 

3.28970 

00 

4.23160-01 

4,00600 

00 

1.00620  00 

-2.22040 

16 

44 

3.15940-01 

3.26930 

00 

4,14780-01 

3,53470 

00 

0,87510-01 

-8.32670 

17 

45 

3.45490-01 

3.24790 

00 

4. 06570-01 

3,02900 

00 

7.60060-01 

-2.22040 

16 

46 

3.75660-01 

3.22580 

00 

3.98520-01 

2.48970 

00 

6.24040-01 

-9,71450 

17 

47 

4.06310-01 

3.20310 

00 

3.90630-01 

1.91790 

00 

4.79690-01 

-1 ,00410 

16 

48 

4,37330-01 

3,17980 

00 

3.82890-01 

1.31470 

00 

3,27330-01 

-1,94290 

16 

49 

4,68600-01 

3,15610 

00 

3,75310-01 

6.81410- 

01 

1,67300-01 

1.38780 

17 

So 

5,00000-01 

3,13210 

00 

3.67880-01 

1.07220- 

01 

1.42220-15 

-6.96630 

19 

29 


NEW  TEST  CASE  1 OF  APERIODIC  TRIDOS 


I 

A(l) 

B(I) 

C(I) 

0(1) 

XHOLD(I) 

OELTA(l) 

1 

0.00000  00 

2.55120 

00 

9.60200-01 

2.91290- 

01 

6.40590-02 

-1,30700 

17 

2 

4,37330-01 

2.60190 

00 

9.60790-01 

5.50590- 

01 

1.30450-01 

-2,77560 

17 

3 

4.06310-01 

2,65190 

00 

9.41760-01 

8,34370- 

01 

1,96970-01 

5,55110 

17 

4 

3.75660-01 

2,70120 

00 

9.23120-01 

1,11770 

00 

2.69400-01 

-4.16330 

17 

5 

3,45490-01 

2,74970 

00 

9.04040-01 

1.40770 

00 

3.41520-01 

-1,24900 

16 

6 

3.15940-01 

2.79710 

00 

8.06920-01 

1.70330 

00 

4,15060-01 

-2.77560 

17 

7 

2,07110-01 

2,64350 

00 

8,69360-01 

2.00330 

00 

4.89760-01 

-6,32670 

17 

a 

2.59120-01 

2.60670 

00 

0.52140-01 

2,30670 

00 

5.65340-01 

-9,71450 

17 

9 

2,32090-01 

2.93260 

00 

0,35270-01 

2.61220 

00 

6.41500-01 

-1.24900 

16 

10 

2,06110-01 

2,97520 

00 

0,10730-01 

2.918S0 

00 

7,17920-01 

-6.93090 

17 

11 

1,01290-01 

3,01620 

00 

0.02520-01 

3.22420 

00 

7,94200-01 

-2,06170 

16 

12 

1,57730-01 

3.05560 

00 

7.06630-01 

3.52010 

00 

6.70230-01 

-1.30760 

16 

13 

1,35520-01 

3,09340 

00 

7,71050-01 

3.02060 

00 

9.45420-01 

*1,66530 

16 

14 

1.14740-01 

3,12950 

00 

7,55780-01 

4.12430 

00 

1.01950  00 

-4,44090 

16 

15 

9.54920-02 

3.16370 

00 

7.40620-01 

4.41370 

00 

1.09210  00 

-2.22040 

16 

16 

7,76360-02 

3,19600 

00 

7.26150-01 

4.69520 

00 

1.16270  00 

-2.22040 

16 

17 

6,16470-02 

3,22640 

00 

7,11770-01 

4,96720 

00 

1.23120  00 

-2.22040 

16 

la 

4,75060-02 

3,25470 

00 

6,97600-01 

5.22030 

00 

1.29690  00 

-2.22040 

16 

19 

3,51120-02 

3.28100 

00 

6.83860-01 

5,47670 

00 

1.35960  00 

-2.22040 

16 

20 

2,44720-02 

3.30520 

00 

6.70320-01 

5,71090 

00 

1,41660  00 

-2.22040 

16 

21 

1.57000-02 

3.32720 

00 

6,57050-01 

5,92930 

00 

1,47410  00 

0.00000 

00 

22 

0.65640-03 

3.34710 

00 

6,44040-01 

6.13020 

00 

1.52520  00 

-2.22040 

16 

23 

3.94260-03 

3,36400 

00 

6.31200-01 

6.31220 

00 

1.57160  00 

-2.22040 

16 

24 

9,66640-04 

3.30020 

00 

6.18700-01 

6,67370 

00 

1,61290  00 

-2.22040 

16 

25 

0,00000  00 

3.39350 

00 

6.06530-01 

6.61310 

00 

1,64070  00 

-2.22040 

16 

26 

9,66640-04 

3.40450 

00 

5.94520-01 

6.72090 

00 

1.67070  00 

-2.22040 

16 

27 

3.94260-03 

3.41330 

00 

5,82750-01 

6.01900 

00 

1.70250  00 

0,00000 

00 

2a 

a. 85640-03 

3,41990 

00 

5,71210-01 

6,00440 

00 

1.71970  00 

-2,22040 

16 

29 

1.57060-02 

3.42440 

00 

5.59900-01 

6.92120 

no 

1,72990  00 

-2.2204b 

16 

30 

2.44720-02 

3.42670 

00 

5.40010-01 

6.92920 

00 

1.73290  00 

-2.22040 

16 

31 

3.51120-02 

3,42690 

00 

5,37940-01 

6.90700 

00 

1,72040  00 

0.00000 

00 

32 

4.75660-02 

3.42510 

00 

5,27290-01 

6.85370 

00 

1,71600  00 

-4.44090 

16 

S3 

6,16470-02 

3,42130 

00 

5.16050-01 

6.76020 

00 

1,69550  00 

-2.22040 

16 

34 

7.78360-02 

3,41550 

00 

5.06620-01 

6.64970 

00 

1,66660  00 

-2.22040 

16 

35 

9.54920-02 

3.40790 

00 

4,96590-01 

6,49730 

00 

1.62920  00 

-2.22040 

16 

36 

1.14740-01 

3.39650 

00 

4.86750-01 

6,31040 

00 

1.50300  00 

-2.22040 

16 

37 

1.35520-01 

3.30740 

00 

4.77110-01 

6.08840 

00 

1,52790  00 

-2.22040 

16 

3a 

1,57730-01 

3,37460 

00 

4.67670-01 

5.83090 

00 

1.46300  00 

-2.22040 

16 

39 

1.81290-01 

3.36030 

00 

4.58410-01 

5.53760 

00 

1,39050  00 

-2,22040 

16 

40 

2,06110-01 

3.34460 

00 

4,49330-01 

5.20640 

00 

1.30810  00 

-2.22040 

16 

41 

2.32090-01 

3.32750 

00 

4.40430-01 

4.64330 

00 

1,21660  00 

-2.22040 

16 

42 

2,59120-01 

3.30920 

00 

4.31710-01 

4.44240 

00 

1.11590  00 

-2,22040 

16 

43 

2.07110-01 

3.20970 

00 

4.23160-01 

4.00600 

00 

1.00620  00 

-2.22040 

16 

44 

3.15940-01 

3.26930 

00 

4. 14700-01 

3.53470 

00 

8,87510-01 

-5,55110 

17 

45 

3.45490-01 

3.24790 

00 

4.06570-01 

3.02900 

00 

7.60060-01 

-9,71450 

17 

46 

3.75660-01 

3,22580 

00 

3.90520-01 

2.40970 

00 

6.24040-01 

*8.32670 

17 

47 

4,06310-01 

3.20310 

00 

3,90630-01 

1,91790 

00 

4^79690-01 

-4.16330 

17 

4a 

4,37330-01 

3,17980 

00 

3.02090-01 

1.31470 

00 

3.27330-01 

-1,66530 

16 

49 

4,66600-01 

3.15610 

00 

3.75310-01 

6.01410- 

01 

i. 67300-01 

1,30780 

17 

50 

5,00000-01 

3,13210 

00 

0.00000  00 

0,36510- 

02 

1.42220-15 

-7,96890 

10 

30 


NEH  TEST  CASE  1 OF  APERIODIC  TRIODV 


I 

A(I) 

8(1) 

C(I) 

0(1) 

XHOLOd) 

DELTA(I) 

1 

0,00000  00 

2.55120 

00 

9.80200*01 

2.91290 

-01 

6,40590*02 

0,00000 

00 

2 

4.3733D-01 

2.60190 

00 

9.60790*01 

5,58590 

•01 

1,30450-01 

-1,38780 

17 

3 

4, 06310*01 

2.65190 

00 

9,41760*01 

8.34370 

•01 

1.98970*01 

1.38780 

17 

a 

3.75660-01 

2.70120 

00 

9,23120*01 

l,il770 

00 

2.69400*01 

*8,32670 

17 

5 

3.45490*01 

2,74970 

00 

9.04840*01 

1,40770 

00 

3.41520*01 

-6.93890 

17 

6 

3.15940*01 

2,79710 

00 

8.86920*01 

1.70330 

00 

4.15060*01 

*9,71450 

17 

7 

2.87110*01 

2,04350 

00 

8.69360*01 

2.00330 

00 

4,89760*01 

*2.22040 

16 

a 

2.59120*01 

2.88870 

00 

8.52140*01 

2.30670 

00 

5.65340*01 

*1,52660 

16 

9 

2.32090*01 

2.93260 

00 

8,35270*01 

2.61220 

00 

6.41500*01 

•2.35920 

16 

10 

2.06110*01 

2.97520 

00 

8.18730*01 

2.91850 

00 

7,17920*01 

*6.43890 

17 

11 

1.81290*01 

3.01620 

00 

8.02520*01 

3.22420 

00 

7.94280*01 

*2.63680 

16 

12 

1.57730*01 

3.05560 

00 

7.86630*01 

3.52810 

00 

8.70230*01 

*1.24900 

16 

13 

1.35520*01 

3.09340 

00 

7,71050*01 

3.82860 

00 

9,45420*01 

9,71450 

17 

14 

1.14740*01 

3.12950 

00 

7.55780*01 

4.12430 

00 

1.01950  00 

*2.22040 

16 

15 

9.54920*02 

3,16370 

00 

7,40820*01 

4,41370 

00 

1,09210  00 

*2.22040 

16 

16 

7.78360*02 

3.19600 

00 

7.26150*01 

4.69520 

00 

1.16270  00 

*2.22040 

16 

17 

6.18470*02 

3.22640 

00 

7.11770*01 

4,96720 

00 

1,23120  00 

-4.44090 

16 

IS 

4.75860*02 

3.25470 

00 

6,97680*01 

5.22830 

00 

1.29690  00 

*2.22040 

16 

19 

3.51120*02 

3.28100 

00 

6.83860*01 

5.47670 

00 

1.35960  00 

0.00000 

00 

20 

2.44720*02 

3.30520 

00 

6.70320*01 

5.71090 

00 

1,41880  00 

•2.22040 

16 

21 

1.57080*02 

3.32720 

00 

6.57050*01 

5,92930 

00 

1.47410  00 

*2.22040 

16 

22 

8.85640*03 

3.34710 

00 

6.44040*01 

6.13020 

00 

1.52520  00 

*2.22040 

16 

23 

3.94260*03 

3.36480 

00 

6.31200*01 

6.31220 

00 

1,57160  00 

*2.22040 

16 

29 

0.86640*04 

3.38020 

00 

6.18780*01 

6,47370 

00 

1.61290  00 

*2.22040 

16 

25 

0.00000  00 

3.39350 

00 

6.06530*01 

6.61310 

00 

1.64870  00 

*2.22040 

16 

26 

9.86640*04 

3.40450 

00 

5.94520*01 

6,72890 

00 

1,67870  00 

*2.22040 

16 

27 

3.94260*03 

3.41330 

00 

5.82750*01 

6.01980 

00 

1.70250  00 

*2.22040 

16 

28 

8.85640*03 

3,41990 

00 

5,71210*01 

6,08440 

00 

1.71970  00 

*2,22040 

16 

29 

1.57080*02 

3.42440 

00 

5,59900*01 

6.92120 

00 

1,72990  00 

-2.22040 

16 

30 

2.44720*02 

3.42670 

00 

5.48810*01 

6.92920 

00 

1.73290  00 

-2.22040 

16 

31 

3.51120*02 

3.42690 

00 

5.37940*01 

6,90700 

00 

1.72840  00 

*2.22(740 

16 

32 

4.75860*02 

3.42510 

00 

5.27290*01 

6,85370 

00 

1.71600  00 

*2.22040 

16 

33 

6.10470*02 

3.42130 

00 

5,16850*01 

6.76820 

00 

1,69550  00 

*4.44090 

16 

34 

7.78360*02 

3.41550 

00 

5.06620*01 

6,64970 

00 

1,66660  00 

*2.22040 

16 

35 

9.54920*02 

3.40790 

00 

4.96590*01 

6,49730 

00 

1.62920  00 

*2.22040 

16 

36 

1.14740*01 

3.39850 

00 

4,86750*01 

6,31040 

00 

1,58300  00 

*2.22040 

16 

37 

1.35520*01 

3.38740 

00 

4,77110*01 

6.08840 

00 

1,52790  00 

•4.44090 

16 

38 

1.57730*01 

3.37460 

00 

4.67670*01 

5.83090 

00 

1,46380  00 

*2.22040 

16 

39 

1.01290*01 

3.36030 

00 

4,58410*01 

5,53760 

00 

1.39050  00 

*2.22040 

16 

40 

2.06110*01 

3.34460 

00 

4.49330*01 

5.20040 

00 

1.30810  00 

*2.22040 

16 

41 

2.32090*01 

3,32750 

00 

4.40430*01 

4.84330 

00 

1.21660  00 

*2.22040 

16 

42 

2.59120*01 

3.30920 

00 

4,31710*01 

4.44240 

00 

1,11590  00 

*2.22040 

16 

43 

2.87110*01 

3.28970 

00 

4,23160*01 

4.00600 

00 

1.00620  00 

*2.22040 

16 

44 

3.15940*01 

3.26930 

00 

4,14780*01 

3.53470 

00 

8.87510*01 

*1 .52660 

16 

45 

3.45490*01 

3.24790 

00 

4,06570*01 

3,02900 

00 

7.60060*01 

*2.00170 

16 

46 

3.75660*01 

3.22580 

00 

3.98520*01 

2.48970 

00 

6.24040*01 

*9,71450 

17 

47 

4.06310*01 

3.20310 

00 

3,90630*01 

1.91790 

00 

4,79690*01 

*1,00410 

16 

48 

4.37330*01 

3,17980 

00 

3,82890*01 

1.31470 

00 

3.27330*01 

*1.94290 

16 

49 

4.68600*01 

3,15610 

00 

3.75310*01 

6.81410 

•01 

1.67300*01 

2.77560 

17 

50 

5.00000*01 

3.13210 

00 

0,00000  00 

8.36510 

-02 

1.42220*15 

*8.72750 

18 

31 


NEM  TEST  CASE  1 OF  APERIODIC  TRIOOM 


TOTAL  DTIME  a 0.154640-01  SECONDS 


1 

A(I) 

8(1) 

C(I) 

0(1) 

XHQLO(l) 

OELTA(I) 

1 

0.00000  00 

2.55120 

00 

9.80200-01 

2.91290 

•01 

6.40590-02 

-1,38780 

17 

2 

4,37330-01 

2.60190 

00 

9,60790-01 

5,58590-01 

1.30450-01 

-1.38780 

17 

3 

4,06310-01 

2.65190 

00 

9.41760-01 

8.34370-01 

1.98970-01 

1,38780 

17 

4 

3,75660-01 

2.70120 

00 

9.23120-01 

1.11770 

00 

2.69400-01 

-6.93890 

17 

S 

3.45490-01 

2,74970 

00 

9.04640-01 

1,40770 

00 

3.41520-01 

-9,71450 

17 

6 

3.15940-01 

2.79710 

00 

8,86920-01 

1.70330 

00 

4.15060-01 

-9,71450 

17 

7 

2,87110-01 

2.84350 

00 

8.69360-01 

2.00330 

00 

4,89760-01 

-2.22040 

16 

S 

2.59120-01 

2,86870 

00 

8.52140-01 

2.30670 

00 

5.65340-01 

-1.52660 

16 

9 

2.32090-01 

2,93260 

00 

8.35270-01 

2,61220 

00 

6,41500-01 

-2,49800 

16 

10 

2,06110-01 

2,97520 

00 

8.18730-01 

2,91850 

00 

7,17920-01 

-6,93890 

17 

11 

1.81290-01 

3,01620 

00 

8.02520-01 

3.22420 

00 

7,94280-01 

-2.63680 

16 

12 

1,57730-01 

3,05560 

00 

7.86630-01 

3.52810 

00 

8,70230-01 

-1.24900 

16 

13 

1.35520-01 

3,09340 

00 

7.71050-01 

3,82860 

00 

9,45420-01 

6.93890 

17 

14 

1.14740-01 

3.12950 

00 

7.55780-01 

4.12430 

00 

1.01950  00 

-2.22040 

16 

15 

9.54920-02 

3.16370 

00 

7.40820-01 

4.41370 

00 

1.09210  00 

-2.22040 

16 

16 

7,78360-02 

3.19600 

00 

7,26150-01 

4,69520 

00 

1.16270  00 

-2,22040 

16 

17 

6.18470-02 

3.22640 

00 

7,11770-01 

4,96720 

00 

1,23120  00 

-2.22040 

16 

16 

4,75860-02 

3,25470 

00 

6.97680-01 

5,22830 

00 

1.29690  00 

-2.22040 

16 

19 

3.51120-02 

3.28100 

00 

6.83860-01 

5,47670 

00 

1.35960  00 

0.00000 

00 

20 

2.44720-02 

3.30520 

00 

6.70320-01 

5,71090 

00 

1,41880  00 

-2.22040 

16 

21 

1,57080-02 

3.32720 

00 

6,57050-01 

5.92930 

00 

1.47410  00 

-2.22040 

16 

22 

8.05640-03 

3,34710 

00 

6.44040-01 

6.13020 

00 

1.52520  00 

-2.22040 

16 

23 

3,94260-03 

3.36480 

00 

6.31280-01 

6.31220 

00 

1.57160  00 

-2.22040 

16 

24 

9.86640-04 

3,38020 

00 

6.10780-01 

6.47370 

00 

1,61290  00 

•2.22040 

16 

25 

0.00000  00 

3.39350 

00 

6.06530-01 

6,61310 

00 

1.64870  00 

0.00000 

00 

26 

9.86640-04 

3.40450 

00 

5,94520-01 

6.72890 

00 

1,67870  00 

-2.22040 

16 

27 

3,94260-03 

3.41330 

00 

5.82750-01 

6,81980 

00 

1.70250  00 

-2,22040 

16 

28 

8,85640-03 

3.41990 

00 

5,71210-01 

6.88440 

00 

1.71970  00 

-2.22040 

16 

29 

1.57080-02 

3,42440 

00 

5,59900-01 

6.92120 

00 

1,72990  00 

0.00000 

00 

30 

2.44720-02 

3.42670 

00 

5,40810-01 

6.92920 

00 

1.73290  00 

-2.22040 

16 

31 

3,51120-02 

3,42690 

00 

5.37940-01 

6.90700 

00 

1,72840  00 

-2.22040 

16 

32 

4,75860-02 

3.42510 

00 

5.27290-01 

6,85370 

00 

1.71600  00 

-2.22040 

16 

33 

6,18470-02 

3.42130 

00 

5. 16650-01 

6,76820 

00 

1.69550  00 

-4,44090 

16 

34 

7,78360-02 

3.41550 

00 

5.06620-01 

6.64970 

00 

1.66660  00 

-2.22040 

16 

35 

9,54920-02 

3.40790 

00 

4.96590-01 

6.49730 

00 

1.62920  00 

-2,22040 

16 

36 

1,14740-01 

3.39650 

00 

4,66750-01 

6,31040 

00 

1.50300  00 

-2,22040 

16 

37 

1,35520-01 

3.38740 

00 

4,77110-01 

6.08840 

00 

1,52790  00 

-2.22040 

16 

36 

1,57730-01 

3,37460 

00 

4.67670-01 

5.83090 

00 

1.46380  00 

-2.22040 

16 

39 

1,61290-01 

3.36030 

00 

4.56410-01 

5,53760 

00 

1.39050  00 

-2.22040 

16 

60 

2,06110-01 

3,34460 

00 

4. 49330-01 

5,20840 

00 

1.30810  00 

-2.22040 

16 

41 

2,32090-01 

3.32750 

00 

4,40430-01 

'4.84330 

00 

1.21660  00 

-2,22040 

16 

42 

2,59120-01 

3,30920 

00 

4.31710-01 

4,44240 

00 

1.11590  00 

-2,22040 

16 

43 

2,87110-01 

3,28970 

00 

4.23160-01 

4.00600 

00 

1.00620  00 

-2.22040 

16 

44 

3,15940-01 

3.26930 

00 

4.14780-01 

3.53470 

00 

8,87510-01 

-8,32670 

17 

45 

3,45490-01 

3.24790 

00 

4.06570-01 

3,02900 

00 

7.60060-01 

-2.22040 

16 

46 

3,75660-01 

3.22580 

00 

3.98520-01 

2.48970 

00 

6.24040-01 

-9,71450 

17 

47 

4,06310-01 

3,20310 

00 

3. 90630-01 

1,91790 

00 

4,79690-01 

-1.80410 

16 

48 

4.37330-01 

3.17980 

00 

3,82890-01 

1,31470 

00 

3.27330-01 

-1,94290 

16 

.4,68600-01 
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